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Chapter I 
INTRODUCTION 


The control of distributed-parameter systems represents 
a real challenge to the system designer, both from a theor- 
etical and practical point of view. Distributed-parameter 
systems arise in various areas. In recent years, due to the 
interest in larger and more flexible structures, the analy- 
sis and design of such systems have steadily increased in 
importance [1-12]. Examples of flexible stuctures are high 
buildings, long bridges, large ships and submersible ves- 
sels, aircraft, rockets, missiles, satellites with flexible 
appendages, astronomical telescopes with large antenas, to 
mention just a few. As a rough approximation, these distri- 
buted-parameter systems, or at least their most essential 
parts, can be considered as vibrating strings, beams, mem- 
branes, plates, shells, and various combinations of these 
components. Such large flexible systems of this nature are 
usually described by a hybrid set consisting of ordinary and 
partial differential equations with two or more independent 
variables [7,9]. The determination of natural frequencies 
and dynamic response of engineering structures, requires 
normally a significant simplification of the actual geometry 
and an approximate method of analysis. 
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Classical structural dynamics problems require the der- 
ivation of the differential equations of motion of typical 
structural members and the determination of the solutio 
ns, where the latter are subject to given boundary condi- 
tions. Such closed-form solutions are limited to relatively 
simple geometries. However, these solutions are of great 
value as they provide an understanding of the dynamical be- 
havior of typical components, such as beams and plates. For 
more complex structures, closed-form solutions are not pos- 
sible and one must consider approximate solutions. Approxi- 
mate methods can provide solutions for a wider range of 
problems and are capable of yielding results of acceptable 
accuracy with reasonable economy of computation. 

The most indicated approximate methods in structural 
dynamics are the classical Rayleigh-Ritz and its finite ele- 
ment version. Substructure synthesis can also be regarded 
as a Rayleigh-Ritz method (13]. The classical Rayleigh- 
Ritz method is less versatile than the finite element meth- 
od, as the latter can be used for more complex structural 
problems [ 13 ] . Complex structures generally represent as- 
semblages of distributed parameter components such as 
plates, beams, columns, etc. Whereas many of these compo- 
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nents may have uniformly distributed properties and by them- 
selves may admit closed-form eigensolutions, when they are 
bounded together in a single structure, the component eigen- 
solution does not have much significance. The eigensolution 
of interest now is that of the complete structure and this 
eigensolution cannot be obtained in closed form. Therefore, 
spatial discretization is unavoidable. 

In contrast to the classical Rayleigh-Ritz method, whe- 
re the expressions for displacements are applicable to the 
complete domain, finite element expressions are defined only 
over a part of the domain, namely an element. It is typical 
of the finite element method, that the resulting discrete 
models possess a large number of degrees of freedom. More- 
over, as with any approximate method, the associated alge- 
braic eigenvalue problem is not an accurate representation 
of the actual system eigenvalue problem, as only the lower 
modes can be computed with satisfactory accuracy and the 
higher modes can be grossly in error. 

In general, the fundamental problem of active control 
of flexible systems is to control a large-dimensional system 
with a much smaller-order controller. Quite often a large 
number of elastic modes may be needed to model the behavior 



4 


of a large flexible spacecraft. However, active control of 
all these modes is not likely to be feasible due to computa- 
tion limitation and errors that arise in modeling the higher 
modes by the discretization process. Clearly, the controls 
must be restricted only to the most significant modes [4]. 
Furthermore, the bandwidths of actuators and sensors prevent 
response to the higher-frequency modes, and these higher 
modes cannot respond to the actual controls. The real just- 
ification for controlling only the lower modes lies in the 
fact that the higher frequency modes in an actual system are 
harder to excite. A conclusion can be drawn at this stage 
that a control system must be designed so as to control the 
modes of most significance, namely the lower modes. To 
guarantee that the controlled modes are sufficiently accu- 
rate, it is often necessary to consider a relatively large- 
order discrete system model. 

The increased order of models representing modern 
structures and the associated control problems, have provid- 
ed the motivation for developing methods permitting analysis 
of such systems, and providing solutions to the control 
problems with little computational effort. Such a method is 
the 'Independent Modal-Space Control' (IMSC) method [6-12]. 
This method is capable of carrying out the control task ef- 



5 


ficiently. Briefly, this method is based on a transforma- 
tion of the system equations of motion to modal space, 
yielding internally independent modal equations of motion. 
Then, the control laws aid designed in the modal space, so 
as to permit independent control • of each individual mode, 
thus providing complete decoupling of the equations of mo- 
tion. This approach allows complete flexibility as to which 
modes to control. The corresponding forces are not actual 
forces but more abstract modal control forces. For control 
implementation, the actual control forces are synthesized 
from the modal forces by a linear transformation. This 
method is easy to implement, especially for high-order sys- 
tems, such as those arising from modeling complex distribut- 
ed-parameter systems. 

Most actuators placement concepts discussed so far in 
conjunction with the IMSC method [14,15], are supposed to 
minimize the part of energy that goes into the uncontrolled 
modes. Indeed, it is shown in Ref. 14 that the work done to 
control the controlled modes does not depend on the actuator 
locations. Of course, this statement holds true if the mode 
participation matrix is nonsingular . This is guaranteed to 
be the case for one-dimensional domains if IMSC is used, but 
cannot be taken for granted in the case of two-dimensional 
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domains, so that care must be exercised in choosing the lo- 
cation of the actuators. 

Another argument can be found in Ref. 15, where a dis- 
cussion on actuators placement is presented. It is remarked 
that actuator locations can be chosen based upon any desired 
criterion which makes the actuators placement as an indepen- 
dent design step. It is also remarked in this reference 
that improper choice of actuator locations can be disastrous 
since the actuator force, or torque, can grow without bounds 
as the actuator locations approch a position for which the 
system becomes uncontrollable. Then, the reference suggests 
a method for optimizing the actuator locations in order to 
minimize the actual control effort, ignoring the problem of 
controllability. In view of this, it is clear that, before 
discussing any minimization aspects, one must place the ac- 
tuators so that the mode participation matrix is nonsingu- 
lar. It should be pointed out that this problem tends to 
disappear as the number of actuators increases. 

The minimum-fuel problem is a very important one in the 
design of space vehicles and various space structures, espe- 
cially for those with lengthy missions. The amount of fuel 
or energy alloted for control is limited to such a degree 
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that it becomes necessary to treat fuel economy as one of 
the predominant factors. In such cases, it is natural to de- 
sign the control system so that it consumes a minimum amount 
of fuel [ 16-23 ] . 

The total fuel consumed during control is measured by 
the time integral of the absolute value of the control vari- 
able. Due to some technical restrictions on the actuators, 
there exist limits on the magnitude of the control forces. 
As an example, the thrust produced by a gas jet actuator is 
limited in magnitude by the saturation of the power ele- 
ments. Controls subject to constraints and minimizing a 
particular functional, such as the consumed fuel, the re- 
quired energy, etc. , are optimal. A control system minim- 
izing the amount of consumed fuel, is referred to as a "mi- 
nimum-fuel system". 

About two decades have passed since the mathematical 
theory of control has been based on powerful variational 
techniques, such as the minimum principle of Pontryagin. 
This method can be used to determine necessary and suffi- 
cient conditions for the control to be optimal [16,17]. The 
necessary condition for optimal control usually specifies 
the nature of the control and the general structure of the 


control system. 
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The derivation of the optimal feedback control law, 
i.e., the explicit dependence of the control on the state 
variables, is a very complex task and in most practical cas- 
es it remains an unsolved problem [18,19]. In these cases, a 
trial and error process is unavoidable for the state deter- 
mination. One may consider such a process as reasonable for 
a low-order system (fourth-order at most), but certainly not 
for a high-crder system, where the computational difficulty, 
is insurmountable. Unfotunately, an accurate model of a 
flexible structure requires a large number of degrees of 
freedom. 

Formulation of the control law in terms of the costate 
variables can be accomplished regardless of the system ord- 
er, but determination of the initial conditions of the cos- 
tate variables for a high-order system is difficult if not 
impossible. The control task is made considerably simpler 
by using the IMSC method where the complexity inherent in a 
high-order system is eliminated. Using this method, the in- 
itial costate variables must be determined for a set of in- 
dependent second-order ordinary differential equations, so 
that treatment of coupled high-order systems is avoided. 
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It is the goal of this study to apply control theory to 
a large flexible structure using the 1 Independent Modal- 
Space Control* (IMSC) method, particularly for minimum- fuel 
control. For the purpose of demonstration, the analysis 
will be focused on a large-order flexible system in the form 
of a plate-like framework. The structure is a beam lattice 
assemblage, as can be seen in Fig. 1.1a. In addition, to 
permit simulation of a gravity-free environment in labora- 
tory, a beam lattice supported by cables, as shown in 
Fig. 1.1b, is analyzed. 

Chapter II contains the problem formulation. In the 
first section, the distributed-parameter system is discre- 
tized via the finite element method for both cases shown in 
Fig. 1.1. Then, in the second section, the eigenvalue prob- 
lem is formulated and modal space decomposition is applied 
to the coupled equations of motion in order to transform the 
coupled set into a decoupled set of differential equations. 
The third section discusses the advantage of using a re- 
duced-order controller, emphasizing the disadvantages in the 
control of the higher modes. In this section, the fundamen- 
tal principles of the IMSC method are introduced, and the 
idea of complete decoupling of the equations of motion is 


established. 





Beam Lattice 

(a) Free in Space 

(b) Cable Supported 
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Chapter III provides a discussion of linear optimal 
control in conjunction with quadratic performance index. 
The object is to suppress rigid body modes and the dominant 
elastic modes, thus providing attitude and shape control of 
a structure in a single control level policy. For control 
implementation, force and torque actuators, such as thrus- 
ters and momentum wheels, are used. 

Chapter IV analyzes aspects of actuators placement and 
their effect on the condition of the mode participation ma- 
trix. The system in this discussion possesses symetric mass 
and stiffness properties in addition to symmetric boundary 
conditions. A proposition stating a sufficient condition, 
related to the singularity of the mode participation matrix, 
is introduced. A detailed proof of the proposition via La- 
place expansion for determinants, concludes this chapter. 

Chapter V analyses the problem of minimum- fuel in 
high-order system. A problem of insurmountable computation- 
al difficulty is made considerably simpler by using the IMSC 
method, where the complexity inherent in a high-order system 


is eliminated. 
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Chapter VI displays numerical examples based on the 
analysis of the beam lattice structure presented in 
Fig. 1.1. 

Chapter VII summarizes the conclusions and provides re- 
commendations for future research. 



Chapter II 
PROBLEM FORMULATION 


2 . 1 SYSTEM DISCRETIZATION VIA FINITE ELEMENT METHOD 


The equations of motion for a flexible structure can be 
derived by means of Hamilton's principle [24]. This deriva- 
tion will lead in general, to a set of hybrid (partial and 
ordinary) differential equations that are difficult to han- 
dle mathematically even for simple structures and it is 
customary to resort to an approximate solution. 

The finite element method is an approximate method of 
analysis which can be used to solve complex structural prob- 
lems. The method consists in taking the displacement mea- 
sures at discrete points in the domain as the unknowns, and 
defining the displacement field in terms of these discrete 
variables. Once the discrete displacements are known, the 
system motion can be described at any given point in the do- 
main by interpolation. In order to achieve system discreti- 
zation, the stiffness and mass matrices as well as the gen- 
eralized force vector for the entire structure, must be 
obtained. In the finite element method we consider the 
structure to be divided into volume elements having finite 
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dimensions and we select certain points in the interior and 
boundary surfaces. The volume elements are referred as "fi- 
nite" elements, since their dimensions are finite, and the 
boundary points are called nodal points or nodes. The vari- 
ous steps in a solution to any problem are: 

1. Discretization of the body by selection of elements in- 
terconnected at certain nodal points. 

2. Evaluation of the element stiffness, mass and force ma- 
trices. 

3. Assemblage of the stiffness mass and force matrices for 
the system of elements-nodes and introduction of dis- 
placement boundary conditions. 

4. Solution of the resulting system equations for the res- 
ponse, due to a given set of initial conditions, under a 
certain design of control law. 

The system shown in Fig. 1.1 will be modeled as a reduced 
dimension model, assembled of flexible beams connected to 
each other by means of rigid blocks, as shown in Fig. 2.1. 
Deriving the potential and kinetic enery expressions as 
well as the non-conservative generalized force, are the 
first steps in achieving system discretization. 
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2.1.1 Potential Energy 

We denote a typical element by j and consider first the 
potential energy expression for the beam member elements. 
From the theory of beams, it is well known that the ratio of 
shear strain energy to bending strain enery is proportional 
to (b/L) 2 , where b is the cross sectional height and L is 
the beam length. Thus, for long thin beams, the shear 
strain energy is very small compared to the bending strain 
energy. For short stubby beams, similar to those given in 
the beam lattice structure, the contribution of those shear 
effects, clearly cannot be neglected and for this reason, 
Timoshenko theory of beams will be considered as a mean of 
accounting for the effects of shear in a simple manner [25]. 
The total slope 3w/3x of the centerline that results from 
shear deformation and bending deformation, can be given as 
the sum of two parts in the following way 


aw (x,t) 

j (x , t) + 


( 2 . 1 ) 


where <#>^(x,t) is the rotation of line elements along the 
centerline due to bending only, while &_.(x,t) accordingly, 
gives the shear angle at points along the centerline. As 
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far as shear strain is concerned, it is assumed that the 
shear strain is the same at all points over a given cross 
section of the beam. That is, the angle Bj(x,t), used here- 
tofore for rotation of elements along the centerline, is 
considered to measure the shear angle at all points in the 
cross section of the beam at position x. Such an assumption 
will greatly facilitate computation. In order to retain 
simplicity and still have the actual shear distribution ef- 
fects, a compensation factor k is introduced and the total 
potential energy for the deformation field of the Timoshen- 
ko-Beam can be written as 



where El , kGA and GJ* are the bending stiffness, the shear 
j 3 J 

stiffness and tortional stiffness, respectively. k is the 
shear distribution compensation factor depending on the 
cross sectional shape and Poisson ratio v, and is given for 
a rectangular cross section as 


lOCHv) 
12 + llv 


( 2 . 3 ) 
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To eliminate t p^ from the total potential energy, we first 
derive the Euler-Lagrange equations and then substitute for 
^ in terms of . Utilizing the principle of mini- 
mum total potential energy, we can derive the equilibrium 
equations associated with the bending and torsion of the 
beam. Hence, let us insert the functional 


F(x,v,w’ .-.p.tp’ ,0,6’) 


El. 


2 kGA . 

+ _1 

dw . 

o * 

2 GJ. 

r d6 / 

2 

dx 

l J 


dx V j 

w J 

+ f 

dx 

L J 


with the three functions w , tp. , 0 into the Euler-Lagrange 

3 3 3 

equation 


d_ 

dx 


3F 




q . = w. ,tp . ,9 . 
1 3 3 3 


(2.5) 


Then, evaluation of Eq. (2.5) for the three functions, re- 
sults in the following equations 


w: h [kGA j (w j ’ " 0 


(2.6a) 


to tEI j^j + kGA j (w j “ *j)] a 0 


(2.6b) 


6: S IGJJ0 j 1 - 0 


(2.6c) 
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Considering elements with constant stiffnesses El., GA , and 

3 j 

GJ* , Eqs. (2.6) yield 

j 

dij>. d^w 

dx = dx^ < 2 - 7 ) 


dw. El. d^w. 

dir " " ' iiGAj dir 1 (2 * 8) 

Substituting Eqs. (2.7) and (2.8) into Eq.(2.2) and using the 
relation G=E/2(l+v), the total elastic potential energy for 
the element can be written as 


V. 

J 


1*1 

El. 

'dV 1 

2 ft 2 

. 12+llv “j. 

d^w . 

, * 
2 EJ. 

, 3 

de . 1 
dx 

J 

r 


2 

dx 

l J 

+ 10 

dx J 

4(1 +v) 

- 


dx 


(2.9a) 


where the first, second and third terms of the integrand 
represent the contribution of bending, shear and torsion to 
the total potential energy. The potential energy for a ca- 
ble member element can be written as 


r x, T 

v j - I ^ (w ' )2 dx 

Jx 3-1 


(2.9b) 


where T Q is the cable tension [13]. 
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In the cable supported beam lattice structure, gravity 
field effect plays an integral role in the expression for 
the total potential energy, and should be added to 
Eqs .(2.9). For this task we will select a typical element 
in gravitation field as can be seen in Fig. 2.2. 

The potential energy can be regarded as resulting from the 
axial force acting through the shortening of the element 
projection on the vertical axis. Denoting a length incre- 
ment along the displaced axis by ds and the projection by 
dy, the potential energy can be written in the form 


r L 

P(y ,t) (ds - dy) 

J 0 


( 2 . 10 ) 


where P(y,t) denotes the axial force due to gravity. Assum- 
ing that the shortening of the projection ds-dy is a small 
quantity of second order in magnitude, and retaining only 
quadratic terms, using the binomial expansion, the shorten- 
ing of the projection can be shown to be 

ds - dy s i (w'(y,t)] 2 dy (2.11) 

where the assumption of small motion was made. Introducing 
Eq . (2.11) into Eq.(2.10) and substituting P(y,t) by the ac- 
tual force acting on the element 




Figure 2.2: A Typical Element in Gravitation Field 
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1 

2 


3 (Qjgyj + Gj ) ( w j ) 2 dy 

J yj_i 


(2.12a) 


where m. is the mass per unit length and G. is an appropiate 
force at the bottom end of a typical element due to the 
weight below. It is clear from the discussion that this 
type of potential energy applies only to the vertical mem- 
bers of the structure, while horizontal members serve as 
load contributors only. For this reason we wish to retain 
both Eqs . (2.9a) and (2.12a) separately, while combining them 
will be essential only when integration takes place in ver- 
tical orientation. 


For the cable elements, the potential energy can be 
written as 


v -± 

j 2 


3 (m gy. + G )(w') 2 dy 

yj-i 3 3 3 


(2.12b) 


where m 

c j 

al, and 


is the mass per unit length for the cable materi- 
is obvious, in view of Eq. (2.12a). 
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2.1.2 Kinetic Energy 

The kinetic energy for a typical element j in a beam 
member can be expressed as 


f x l 

m . 

'3w^' 

2 J. 

r 3w f / 

2 j . 

30 . 

2* 

J 



+ -1 

1 

+ CLl 

1 



2 

3t 

0 

4m 

3 t 

' 2 

3t 


Jx j-1 

- 

l J 


l j 


t J 



where m. represents the mass per unit length, J. is the mass 
J J 

moment of inertia of an element of unit length about the 
bending axis and J 0 ^ is the polar mass moment of inertia of 
an element of unit length. Moreover 3w^/3t is recognized as 
the translational velocity, 3wj/3t as the angular velocity 
of the element about the bending axis and 36^/3t as the an- 
gular velocity of the element about the torsion axis. Ac- 
cordingly, the first term of the integrand represents the 
translational the second the rotational and the third the 
torsional kinetic energy. The kinetic energy for the rigid 
blocks, can be written in view of the notations of Fig. 2.1c 
as 


T. 

J 




xG 




e 2 . 

yj 


(2.13b) 


where M is the mass of the rigid block and I and I 

xG yG 

its mass moment of inertia about its lateral axes x and 


are 


Y> 
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respectively. For the cable elements, the kinetic energy 


can be written as 


f y-i m . 

T. = I ^ (w.) 2 dy (2.13c) 

^ ) v ^ 

y j-l 


2.1.3 Nonconservative Generalized Force 

In the derivation of the discretized system equations 
of motion, we are making use of the virtual work concept in 
order to obtain the nonconservative generalized force ex- 
pression 

oW . = P [f , n>! m ] 5 (w w' 9 ] T dx (2.14) 

j J wj wj J 3 3 

x j-l 

where f ., m » , and m Q . are the distributed force and two 
wj wj 0 j 

moment components associated with the displacements w^ , w! , 
and 9^, respectively. 
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2.1.4 Finite Element Model 


See also Meirovich [13], chapter 9. 


Referring to Fig. 2.3a we express the displacement com- 
ponents at any point inside the element j in the form 


w (5,t) = L^(C) w -i (t) 

J -W ~J 

0j(€.C) = Lg(C) 6j(t) 


where 

L u (?) - (!•„/?) L W3 ( ? ) L Uu ( 5 )] T 


(2 . 15a) 


(2.15b) 


(2.16a) 


p 9 (5) - [L 3i ( 5) L 92 (5)] T 


are vectors of 
Wj(t) = [w^OO 

0j(t) = Ce^._ 1 (t) 


interpolation functions and 
hwj_ L (t) Wj (t) hw^(t)] T 

QjCt)] 1 


(2.16b) 


(2 . 17a) 


(2.17b) 


are vectors of nodal displacements where h is the length of 
the element. To satisfy compatibility, it is necessary that 
the interpolation functions vector be continuous up to the 
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derivative of one order lower than the highest derivative 
appearing in the associated differential equation. Most of 
the element formulations that have been developed are based 
on polynomial expansions. Therefore, for a beam in bending 
represented by a fourth order differential equation, cubic 
polynomials are admissible, while for the torsion part that 
is represented by a second order differential equation, li- 
near polynomials are admissible. For the discretization 
process at hand, we select the following interpolation func- 
tions 

L Ul <0 - H 2 - 25 3 , l V2 (C> - 5 2 - C 3 

V«> ' 1 - 352 + 253 • \ <«)-'« + 2 5 2 - 5 3 (2 - 18a) 


• L 6 2 (5) ' 1 - 5 


(2.18b) 


For a member with rigid blocks, the following relations can 
be derived in view of Fig. 2.3b. 


w (C,t) - L^(£) w^(t) = L*^(£) w (t) 


(2.19) 


Vi - Vi + Vi f 


W = cj T 

j-i j-i 


(2.20a,b) 


w. 


W J - W j 2 


W = XJ 1 

j j 


(2.20c,d) 
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Referring to Eqs. (2.16a) and (2.17a) and substituting 
Eqs . (2.20) into (2.19), we obtain 

w.(t) = C R w.CO (2.21) 

where C.^ is the rigid constraints transformation matrix 



( 2 . 22 ) 


Substituting Eq. (2.21) into Eq.(2.19), the modified interpo- 
lation function vector L*(£) becomes 

- W 

LJ(0 = cj L y (0 (2.23) 

We will use two reference systems, local and global. For a 
single element, nodal variables refer to local numbering 
system, but for an assembly of elements we shall refer to 
global system. The relation between global and local coor- 
dinates x and £ respectively, can be concluded from Fig. 2.3 



dx = -h d£ 


(2.24) 
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The element beam matrices can be obtained in terms of the 
local coordinate £ by substituting Eg. (2.24) into (2.9a). 
The potential energy for a typical element j in the beam 


lattice has the form 


v! 


rl l 

El — 

'd 2 w (C ,t) 

2 

, 12+llv 

fil 

~d 3 v.a,c) 

Jo 

^ j h 4 

<£ 2 

+ 10 

A. h 6 
J 



EJ* 

2 

"de 1 (? ,c)‘ 

T 

4(l+v) j 


, d 




(2.25a) 


Inserting Eqs.(2.15) into (2.25a), one obtains 


v j -H<'> \ } ?j<'> + T§>> Sj 9l (t) 

where 

l d2 i w <«> < 2 £cO 

■ h3 j 0 EI j -7i W~ « 

. i r 1 12+iov EI * d V?> d3 ^< E) 

h* J 0 10 Aj — dp dp - d5 



1 

h 


r l EJ* dL e (0 dL*(0 
, 0 4(l+v) ~ d l d ~ 


(2.25b) 


(2.26a) 


(2.26b) 


are element stiffness matrices associated with the tran- 
sverse and torsional motions, respectively. Using 
Eqs.(2.16) and (2.18), the element stiffness matrices become 
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‘ 12 6 -12 6 


" 144 72 -144 72 ' 

El 

4-6 2 

12+llv El 2 

36 - 72 36 

= h? 

12 -6 

10 Ah 3 " 

144 -72 


syram. 


symm. 


4 


36 


(2.27a) 


K 


9j 


EJ* 

4 (l+v)h 


1 

-1 


-1 

1 


(2.27b) 


where all parameters are assumed to be constant throughout 
the element. The potential energy for an element j in a ca- 
ble suspension is obtained by substituting Eg. (2.24) into 
(2.9b). The result is 


V. 

J 


f r i, 

dWj (C,t)' 

'o 

dC 


dC 


(2.28a) 


Using Eg. (2.15) 


V. 

J 


= J K c . w.(t) 


where 



.^f 1 

h. |„ 


T 

c 


d_L Q (0 dL^(0 
~ d£ ~ 


dt. 


(2.28b) 


(2.29) 
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using the linear interpolation functions Eqs. (2.16b) and 
(2.18b), the cable element stiffness matrix becomes 


K . 
cj 



-1 


1 


(2.30) 


Following this pattern, the additional stiffness matrix due 
to gravity contribution can be obtained from Eq. (2.12a) 


V. 

J 


- f 1 

2 Jo 


(m^g (I- p ) h + ] 


,t) 


dC 


d? » j wT(t) K w. (t) 
L '] * J ~j 


(2.31) 


where 





dL (?) dL^(C) 

- W 

d? d? 


dC 


|»1 

®.g5 

J 0 J 


dL (5) dljtt) 

W W 

d? d5 


dC 


(2.32) 


Using Eqs. (2. 16) and (2.18) 
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36 0 


1 , G 

30 m § h 


symm. 


(2.33) 


The potential energy due to gravity contribution for element 
j in the cable suspension, is similar to Eq. (2.28), or 


v i Vi 


(2.34) 


Since cables in transverse vibration and beam in torsional 
vibration represent entirely analogous dynamical systems 
from a mathematical point of view, the same interpolation 
functions can be used for both cases. Therefore, K cG ^ in 
Eq.(2.34) has the form 


If G ] dL (O dL x (0 

0 cj § h c d£ d£ 


dL U) dL*(0 

» cj * 5 -V (2 ‘ 35) 


Using Eqs . (2.16b) and (2.18b) 


3 . G 

K CGj = 2 m c g + TT 


1 -1 
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A similar pattern is adopted to obtain the element mass 
matrix, where the assumption of a member consisting of one 
element is made. Substituting Eq.(2.24) into (2.13a) 


T. 

2 


h 

2 


[w. (£,t) ] 2 + j. 
J J J 


1 

h 2 


f ~dw(£ ,t)1 ' 

L dZ 1 


+ J 


Oj 




dZ 


(2.37a) 


inserting Eqs,(2.19) and (2.23) into Eq. (2.37a), one obtains 


T. 

J 


I M 9J 8 j( t) 


(2.37b) 


where 


M . 

wj 


1 r 1 T d V» 

hjn i dC 


= h m L , (5) L (Z) dZ + 

n 2 w ~w 


dZ 


dZ 


(2.38a) 




LgCO L T (5) d£ 


(2.38b) 


T 

CM .C_ and M. . are element mass matrices associated with 
R wj R 0j 

the transverse and torsional motion. Using Eqs.(2.16) and 
(2.18), the element mass matrices become 
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156 22 54 -13 ' 


' 36 3 -36 3 ' 


mh 

4 13-3 


4-3-1 

. 

420 

156 -22 

30h 

36 -3 



symm. ^ 


symm. ^ 

J 


(2.39a) 



J h 



1 

2 


(2.39b) 


where all the parameters were assumed to be constant 
throughout the element. The kinetic energy for element j in 
the cable suspension is obtained by substituting Eq.(2.24) 
into (2.13b), so that 


T. 

J 


h r 1 

m (w.(C,t)] 2 d£ 
i o 3 J 


(2.40a) 


Using again the linear interpolation functions Eqs. (2.16b) 
and (2.18b), we obtain the cable element mass metrix 


T. 

3 


T "j (t) M C j ^j (t) 


(2.40b) 


where 
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m h 
c c 


6 


2 

1 


1 

2 


( 


The discrete element force vector can be obtained 
Eq.(2.14). Using local coordinate ? 


6 W. = h 

J 


Jo 


f . OW. (£ ,t) - m >. T- 6 
wj j wj h 


dw^ (? ,t) 


d? 


+ m ej 59^5 




6C 


Substituting Eqs.(2.15) into (2.42a) 


,1 


5W. = h f . l/-(£) C 0 6w. (t) d? - 
J j 0 WJ -w R ~3 


"„'j — — C R <t)d5 


+ h 


m Lt(5) C- 50. (t) d? 
Q 0 j -0 K. J 


- 4 + ?ej 


(2 


where 


F . - h f . Cl L (?)d? 
-wj J Q wj R -w 


_ dL (?) 

T j 

m c: d? 


Vj R d? 


(2 


r 1 

i . L, 

Jo 93 '< 


F 9j - h . " M « 


(J 


.41) 


from 


.42a) 


.42b) 


.43a) 


.43b) 


are the element generalized force vectors. 
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We will regard the forcing functions as concentrated 
forces and moments acting at the element nodes. These forc- 
es and torques can be represented as distributed forces and 

torques by writing them in the form F 6(5-5 ) and M 6(5-5 ) 

J j 

where 6(5-£ ) and 6(5-5 ) are spatial delta functions. Sub- 
i j 

stituting these expressions into 2qs.(2.43), we can write 


F . - [h F. , M . , h F M ] T (2.44a) 

~ w 3 J-l wj-1 n f j wj J 



th Vi hM 63 lT 


( 2 . 44b ) 


where , M w j =1 , Mqj.i and Fj, M w j, M 0 j are concentrated 

forces, bending moments and torques at the j-l and j nodes 
respectively. 


The potential energy, kinetic energy and virtual work 
for member s are obtained by summing up the contributions of 
the individual elements, in a process known as " assembling" . 
Two types of members participate in the assembly, beams and 
cable supports which are regarded as 4th and 2nd order mem- 
bers respectively, according to the degree of their differ- 
ential equations of motion. With this notation it will 
prove convenient to introduce a 3(m 4 +l) and a (m 2 +l) dimen- 
sional displacement vectors for a 4th and a 2nd order mem- 
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ber, where m 4 and m 2 are numbers of elements in 
of members, respectively. For a beam-member 

U = fw n hwA Wl hwi ... W hw' 0 n 01 

-S4 U u 1 1 m 4 m4 u 1 

and for a cable-member 


these type 


m4 


] 


T 


(2.45a) 


u s2 - t w 0 W 1 



(2.45b) 


The member stiffness matrix can be written in a block diag- 
onal form as follows. Starting with a beam member, we have 


K 


S4 


K ! 0 

W I 

4 


K. 


3(m 4 +l) 


(2.46) 


where each block can be schematically displayed as follows 



(2.47a) 


2(m 4 +l) 
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in which the shaded blocks indicate summation of the over- 
laping entries. We notice the distinction between a verti- 
cal member and a horizontal member in the beam lattice. The 
former is affected by gravity, while the latter serves just 
as a load contributor. Therefore, for a vertical beam mem- 
ber, using Eqs.(2.27) and (2.33), we obtain 
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144 

72 

-144 

72 





72 

36 

-72 

36 





-144 

-72 

288 

0 

-144 

72 



72 

36 

0 

72 

-72 

36 


El 2 



-144 

-72 

288 

0 

, 12+llv 






+ 10 

Ah 2 



72 

36 

0 

72 


288 

0 

-144 

72 

0 

72 

-72 

36 

-144 

-72 

144 

-72 

72 

36 

-72 

36 



I 


, G 
m2 + h 


' 36 

3 

-36 

3 



3 

4 

-3 

-1 



-36 

-3 

72 

0 

-36 

3 

3 

-1 

0 

8 

-3 

-1 



-36 

-3 

72 

0 



3 

-1 

0 

8 


0 


+ 


72 

0 

-36 

3 

0 

8 

-3 

-1 

-36 

-3 

36 

-3 

3 

-1 

-3 

4 


0 



40 



36 

0 

-36 

6 



0 

6 

0 

-1 



-36 

0 

72 

-6 

-36 

6 

6 

-1 

-6 

8 

0 

-1 



-36 

0 

72 

-6 



6 

-1 

-6 

8 


0 


0 


(2.48a) 


72 

-6 

-36 

6 

-6 

8 

0 

-1 

-36 

0 

36 

-6 

6 

-1 

-6 

2 


1 -1 

-1 2 -1 


K 


9 


EJ* 

4 ( 1+v ) h 


-1 2 


0 


0 


2 -1 

-1 1 


(2.48b) 


Here, we should note that for a horizontal beam-member, only 
the first two matrices in Eq. (2.48a) are to be included. 
Finally, for a cable-member, the stiffness matrix can be ob- 
tained using Eqs.(2.30) and (2.36), with the result 
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1 -1 


-1 2 -1 


K 

s 2 


3 

2 m C S 


+ 


t c + g1 


h 


-1 2 


0 


0 


2 -1 

-1 1 


(2.49) 


A similar pattern can be hold to obtain the mass matrices 
for a beam-member as well as for a cable-member. 


For the member force vector, we can write 
- F s 4 - [ - F w - 0 ^ 3 ( 04 + 1 ) (2.50) 

where F and F are 2(m 4 +l) and (m 4 +l) vectors respectively, 

~ W ~ 0 

obtained from Egs.(2.44) in a manner similar to the member 
stiffness and mass matrices. 


Making use of the above matrices and vectors, the mem- 
bers equations of motion can be written in the matrix form 


M u + K u = F 
S 4 ~ s 4 s 4 - s 4 " s 4 


S4 a 1)2) •• • 


(2.51a) 


M u + K u _ = 0 

s 2 ~ s 2 s 2 ' s 2 


S2 1 , 2 1 . . ■ 1^2 


(2.51b) 
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for a 4th and 2nd order members respectively. M 4 and M 2 
represent the number of 4th and 2nd order members in the 
structure. It will prove convenient to consider the total 
dimension of the beam lattice according to the number of 
beam members along each side. Let n^ and represent these 
numbers along the lattice sides a and b respectively, and 
n c , the number of cable members in each cable branch. 


" n a (n b +1) + n b (n a +1) (2.52a) 

M 2 - 2n fi (2.52b) 


Equations (2.51) represent a set of M=M 4 +M 2 disjoint equa- 
tions for the independent members in the entire structure. 
For the assemblage of discrete members, the boundary dis- 
placements must match. The displacements at the nodes shared 
by several members must be the same for every such member 
and corresponding forces should be statically equivalent to 
the applied forces. Note that displacements include rota- 
tions and forces include torques. At this point let us in- 
troduce the N^=3M 4 (m 4 +1 ) +M 2 (m 2 +l ) dimensional disjoint dis- 
placement and force vectors 


“d " ["41 "42 ••• “I M “21 ^22 


y 2 


m 2 


(2.53a) 


d 






(2.53b) 
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With these expressions, the equations of motion for the dis- 
joint structure, have the reduced matrix form 


M, a, + K. a , => F, (2.55) 

d-d d ~d ~d 

It is clear that the disjoint vector u d contains local coor- 
dinates some of which are redundant for the description of 
the assembled structure. It will prove convenient to intro- 
duce a new displacement vector q representing the nodal dis- 
placements of the complete structure in terms of components 
along the global system coordinates. The total number of 
nodes in the entire structure can be determined in terms of 
n , , n c , and m 4 ,m 2 ,M 4 ,M 2 

N 4 = (n +l)(n b +l) + M 4 (m 4 -1) , N 2 - 2a, - 2 + M 2 (m 2 -1) (2.56a,b) 

In view of this number, the order of the independent global 
coordinates vector q is 


N = 3[(n +1) (n, +1) + M 4 (m 4 -1) + [2n - 2 - M 2 (o 2 -l)] 

3D C 


(2.57) 
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The relation between the vectors u d and q can be written in 
a matrix form as 

a d (t) = C q(t) (2.58) 

C is a N d xN^ rectangular matrix reflecting the constraints 
expressions, relating the disjoint structure in the local 
coordinate set to the assembled continuous structure in the 
global coordinate set. Introducing Eq.(2.58) into Eq.(2.55) 
and premultiplying by C T , the structure equations of motion 
will take the matrix form 

Mq(t) + K q(t) = Q(t) (2.59) 

where 

M = C T M d C , K = C T K d C , Q - C T F d (2.60a,b,c) 

and 

,T 

q = tw X 9 X1 e yi w 2 0 X2 e y2 ... 0^ 0 yN(+ W N ^ +1 ... w Nu+N2 j 

( 2 . 61 ) 

q is the nodal displacement vector and Q represents the vec- 
tor of nodal control forces. The elements in the nodal dis- 
placement vector, w^, 0 xi , 0 y ^ are the nodal global coordi- 

nates, i.e., lateral displacement and two rotations 
respectively, in the ith node. To add the rigid blocks con- 
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tribution to the kinetic energy, we can associate the veloc- 

• • • 

ity components of Eq. (2.13b) Wj , 8 x j, 8 y j , with the appropri- 

• • • 

ate nodal velocities q , q , q , respectively, and sim- 

i i+1 i+2 

ply add the mass and two moments of inertia of these blocks 
to the related entries in the main diagonal of the system 
mass matrix. 


2.2 ' THE EIGENVALUE PROBLEM AND MODAL SPACE DECOMPOSITION 

The object is to reduce the transverse motion of the 
structure to zero. This can be done most effectively by ac- 
tive control techniques. Customarily, controls are designed 
by working with the state equations of motion. In the 
state- space framework, methods have been developed for state 
and output feedback problems which have been proved to be 
computationally’ -attractive. It turns out that for a discre- 
tized structure, as described by Eq.(2.59), a configuration 
space approach is simpler and more useful. 

To achieve the control task, the IMSC method is em- 
ployed. As a first step, it is necessary to solve the ei- 
genvalue problem associated with Eq.(2.59), where the latter 
can be written as 

<4 M 3 r = K 9 r 


r = 1,2, ...,Nq 


(2.62) 
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M and K are real symmetric and positive definite matrices, 
so that all the eigenvalues are real and positive, where 
the u r ' s represent the natural frequencies of the system and 
q r are the associated orthogonal set of eigenvectors. The 
eigenvectors can be normalized by setting 


T 

q M q 
~r -i 


6 

rs 


T 

Sr 


K S r 


2 - 
U 4 0 

r rs 


r,s = l,2,...»Nq (2.63a,b) 


where 5 rs is the Kronecker delta. Arranging the eigenvec- 
tors in a N xN modal matrix 

q q 

u = [qi q 2 ••• S Nq ] (2.64) 


Eqs.(2.63) can be rewritten in the compact matrix form 
U T MU = I , U T KU = A (2.65) 


where I is the identity matrix of order N and A is the 


q 


diagonal matrix of eigenvalues. 

Using the expansion theorem [24], the response of the system 
at any time can be expressed as a linear combination of the 
eigenvectors. Defining £(t) as a dimensional modal space 
vector, we" introduce the linear transformation 


q(t) = U£(t) 


( 2 . 66 ) 
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Substituting Eq.(2.66) into Eq.(2.59), premultiplying by U T 
and using Eq.(2.65), we obtain the modal equation 

§(t) + A £(t) = N(t) (2.67) 

where 

N(t) = U T Q ( t) (2.68) 

is the Nq - vector of generalized control forces associated 
with the modal vector £(t). Eq.(2.67) represents a set of 

modal equations 

M C) + w r = N r (t) r=l,2,...,Nq (2.69) 

Although Eqs.(2.69) appear as a decoupled set, modal coordi- 
nate recoupling occurs in general through the control forc- 
es, as can be recognized from Eg. (2.68). To achieve a com- 
plete decoupled form, each modal control force N r must 
depend on the rth mode only, which is the essence of the 


IMSC method. 
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2.3 CONTROL OF THE DOMINANT MODES 


The fundamental problem of active control of a flexible 
system, is to control a large-dimensional system with a much 
smaller-order controller. Although a large number of elas- 
tic modes may be needed to model the behavior of a large 
flexible structure, active control of all of these modes is 
out of question due to computation limitations and errors in 
the higher modes because of the discretization process. 
Clearly, the controls must be restricted to the most signi- 
ficant modes, which are the lower modes. 

The justification for controlling lower modes only lies 
in the fact that higher frequency modes are very difficult 
to excite. Furthermore, the bandwidths of actuators and 
sensors cannnot respond to the higher frequency modes. 
Hence, although Eqs.(2.69) permit independent control of all 
the discretized system modes, we will concentrate in cont- 
rolling only a number of dominant modes. This number can be 
determined by examining the modes participation in the res- 
ponse. To accomplish this, the eigensolution will be rear- 
ranged according to the increasing order of the eigenvalues 
and only the first c of Eqs.(2.69) will be retained for the 


control task. 
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Perhaps this is the point to emphasize the unique ad- 
vantage of the Independent Modal Space Control method. Aft- 
er the system equations have been transformed to the modal 
space, the control laws will be designed in the modal space, 
permitting independent • control of each individual mode. 
This approach allows a complete flexibility as to which mode 
to control. For control implementation, the actual control 
forces are synthesized from the modal control forces by a 
linear transformation, as we shall see latter in this sec- 
tion. It appears that this approach is simpler to imple- 
ment, especially for high-order system such as those arising 
from modeling complex distributed-parameter systems. Let us 
assume that only a subset of the modeled modes will be cont- 
rolled and use the subscripts C and R to denote the cont- 
rolled and residual (uncontrolled) modes, respectively. 
Then, introducing the partitioned forms 


§<c> 



» - '»= i V 

(2.70a,b,c,d) 


Equation (2.67) can be separated into the two equations 
L + A c §c = ?c (2.71a) 


§R + A R Sr “ Sr 


(2.71b) 
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where 

- “c 9 . ?R - U R 9 


(2 . 72a, b) 


are the partitioned parts of the control vector. 

Substitution of Eq. (2.70a) into Eq.(2.66), yields the co- 
nfiguration space displacement vector 


q(t) - U c « c Ct) + U R §R (t) 


(2.73) 


Assuming that the contribution of the uncontrolled modal 

vector £ (t) to the actual response is negligible, we can 
R 

truncate Eq.(2.73) by ignoring U £ (t). 

R — R 


For control implementation, thruster and torquer ac- 
tuators will be placed at some nodes, the number of which 
equals to the number of controlled modes. The control forc- 
es are first designed as generalized forces in the modal 

jmssass^^ 

space, and then synthesized to form the actual actuator 
forces . 

The discrete forcing vector Q(t) in Eq.(2.59) will be gener- 
ated by using c discrete actuators 

Q(t) - B F(t) (2.74) 

where F(t) is a c-dimensional vector of thrusters and tor- 
quers, while B is a full rank N^xc modal participation ma- 
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trix. Substituting Eq. (2. 74) into Eq. (2.72a), the expres- 
sion for the actual actuators forces is 

F(t) = (U^B) -1 N c (t) (2.75) 

where the actuators are placed so that (U B) is a nonsingu- 

L 

lar matrix. 

Substituting Eqs.(2.74) and (2.75) into Eq. (2.72b), the re- 
sidual modal force vector is obtained as 

N R (t) = U R B(U* B)" 1 N c (t) (2.76) 

We can summarize now the modal equations of the cont- 
rolled and uncontrolled modes as follows: 

£ (t) + A e (t) » N (t) (2.77a) 

- c c - c -c 

'L(t) + a 5 (t) = u* b(u* B)" 1 n (t) 

- R R -R R c -c (2.77b) 

Equations (2.77) illustrate clearly that the c entries of 
the modal vector £ c (t) can be controlled independently, whi- 
le the R entries of the residual part of the modal vector 
£ R (t) are excited by the spillover of the control forces 

into the uncontrolled modes. 



Chapter III 


LINEAR OPTIMAL CONTROL WITH QUADRATIC 
PERFORMANCE INDEX 


An optimal control system is defined as one in which a 
certain performance index is minimized. Selecting the prop- 
er performance index for a complex control system is in gen- 
eral a difficult task. As stated in the previous chapter, 
independent control of the modes is achieved if the control 
force associated with each modal coordinate, depends only on 
this coordinate. This guarantees complete decoupling of the 
modes. In view of this approach, we can determine a perfor- 
mance criterion for each mode independently. Because every 
modal coordinate is controlled independently of any other 
modal coordinate, the overall system cost will then be the 

sum of all modal costs 
c 


J = 
c 



(3.1) 


Evidently J c is a minimum if every J r in the summation is a 
minimum, because minimization of each term can be carried 
out independently. At this point we abandon the configura- 
tion space expression from chapter II and replace the sec- 
ond-order modal equations Eqs. (2.71a) by first-order modal 
state equations. 

Let us consider the modal state variables 
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This type of state definition is different from the one used 
in the IMSC method [26]. In Ref. 26, the modal velocity is 
eigenvalue-dependent such that the modal matrix A r is skew 
symmetric in « r and vanishes from the modal equations for 
possible rigid body modes. To handle the entire system in 
presence of rigid body modes, dual- level control scheme was 
introduced, in which the first level is a proportional con- 
trol providing control to the rigid body modes. This con- 
trol provides "artificial" stiffness rendering the stiffness 
matrix positive definite. Second level controls 'are then 
designed to provide final control to the complete system. 

The state defined by Eqs.(3.2) permits one control policy 
for the entire system, including control of the rigid-body 


modes . 
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In view of the state Eqs . ( 3 . 3 ) , we consider the follow- 
ing performance indices 



(3.5) 


where t f is the final time, Q r are 2x2 weighting matrices 
and R r are weighting factors. For the rigid body modes, Q r 
are chosen independently of the eigenvalues, or 



r = 1,2, ... ,n 


(3.6) 


For the elastic modes, Q are chosen as 

56 r 



r = n +1 ,n +2 , . . . ,c 
r r 


(3.7) 


In view of Eq.(3.7) the first term of the integrand in 

Eq.(3.5) represents the total energy associates with the 

controlled elastic modes. R is the modal control effort 

r 

expenditure during the interval (0,t ). Increasing R^ , 
places a havier penalty on the control effort. 

The minimization of the modal cost J r leads to the optimal 
controls [16,17] 

N r (t) - - R" 1 B* K r x r (t) r = 1,2, ... ,c (3.8) 
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where K r (t) is the 2x2 symmetric matrix satisfying the modal 
differential matrix Riccati equation 


K = - Q + K B R _1 B T K 
r r r r r r 


- K A - A T K 
r r r r 


(3.9) 


r = 1,2, . . 


The steady-state solution is of great interest and will be 

* 

obtained by setting K r =0 in Eqs.(3.9). This yields the 
algebraic matrix Riccati equation. With the notation 



K 

K 

K = 

llr 

1 2r 

r 

K 

K 


2 lr 

22r 


r = 1,2,. ..,c (3.10) 


the solution to the gain coefficients K ilr , K t 2r / K 22r can 
be shown to be 


"nr" 


W 22r - 2R r “P«Ur + R c “r> + 2R r‘ < 3 ' lla > 


-,4 


K 22r * ! R r «22r ' 2R t “I + 2R r (R r < + \ '.l, 1 ' 


(3.11b) 


K 


12r 


= - R w 2 + (R 2 id 4 + R Q. . ) 
r r r r r x llr 


(3.11c) 


T ly2f«»*yC 
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where Q llr , Q 2 2r are general coefficients of the diagonal 
weighting matrix Q r . Substituting Eqs.(3.11) into Eqs.(3.8) 
and using Eqs.(3.2), the optimal modal controls are obtained 
in the form 


K 


N r (t) 


12r „ 


22r ; 


R 'r 
r 


5 (t) 


R ’r 
r 


r l,2,.».,c 


(3.12) 


Insertion of Eqs.(3.12) into Eqs.(2.71a) results in the sec- 
ond-order differential eouations 


K 22r • 

5 r co + 5 r (t> + 


K 

2 , 12r 


oj“ + 
r 


5 r (t) 


= 0 r = 1,2,. ..,c (3.13) 


The general solution to Eqs.(3.13) is known to be 


? r (°) 

5 r CO X 

r a - 4>* 


5 (0) -G ai c 

e r r cos (a), t - ) + e sinCa^t) 

dr r w, 

dr 


Mr 


where 


» 1,2,. ..,c (3.14) 


K, „ 

9 o , 1 2 r 

“nr * “J + ~R~ 


(3.15a) 


K 


22r 


2(R 2 oj 2 + R K ) : 
r r r 12r 


(3.15b) 


(i), — oj (l - ; 2 ) 2 

dr nr r 


(3 . 15c) 
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V 

= arct an j ( 3 . 1 5 ci ) 

r U -;*) 1 

Eqs.(3.14) imply that l-c*>0, which places a lower bound on 
the control factors R r . Using Eqs . ( 3 . lib, c ) and (3.15b), we 
obtain 


R > 


Q 2 

y 22r 


r 4 (0 + Q_ „ ^T 

ulr x 22r r 


r » 1,2, . .. ,c 


(3.16) 


From Eq.(3.16), in view of Eq.(3.6) and (3.7) 


> 4Q 


22r 


r = 1,2, . . . , n (3.17a) 


llr 


for a rigid body mode and 


R r =• ii? 


r * n r +l ,ri r +2 , . . . , c (3.17b) 


for an elastic mode. 

Substituting Eqs . ( 3 . 11b, c ) into Eq. (3.15b), and setting 
Q llr =w| and Q 2 2 r =1 ' we obtain 


1 

2 


l l 

1 - 2p + 2p (1+p) 2 


p 2 (1 + p) 


r = 1,2, . . . ,c 


(3.18) 


where p is a dummy variable 
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p = R oj^ r * 1,2, ... ,c (3.19) 

v r r 

Figure 3.1 provides a plot of the modal damping factor i r vs 
p. Using this plot and in conjunction with Eq.(3.19), one 
can select the desired control weighting factors R r 
( r=l , 2 , . . . , c ) to meet certain modal response performance 
criteria, such as settling time, etc. A numerical example 
for linear optimal control of the beam lattice is presented 
in Ch. VI. 




Chapter IV 


THE EFFECT OF ACTUATORS PLACEMENT ON THE MODE 
PARTICIPATION MATRIX 


Most actuators placement concepts discussed so far in 
conjunction with the IMSC method [14,15] are directed toward 
minimizing the energy that goes into the uncontrolled modes. 
These concepts provide guidelines for achieving this goal. 
However, there are some other constraints that must be sa- 
tisfied in conjunction with the minimal solution. 

It has been shown that the work done to control the 
controlled modes does not depend on the actuator locations. 
Of course, this statement holds true if the mode participa- 
tion matrix is nonsingular . This is the case for a one-di- 
mensional domain if IMSC is used. This cannot be taken for 
granted in the case of two-dimensional domains, so that care 
must be exercised in choosing the location of the actuators. 
In view of this, it is clear that, before discussing minimi- 
zation aspects, one must place the actuators so that the 
mode participation matrix is nonsingular. 

In this chapter we establish a sufficient condition for 
the mode participation matrix U^B that asserts its singular- 
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ity, where U T is the cxn upper part of the modal matrix and 
C 

B is a full rank nxc actuators placement matrix, in which n 
is the order of the system and c is the number of controlled 
modes. The analysis provides a method of recognizing unde- 
sirable actuators placement configurations. 

We start by examining the structure of the modal matrix 
U. The nxn modal matrix has the modal vectors as its co- 
lumns. Furthermore, the modes are arranged in increasing 
order of magnitude of the associated eigenvalues. We con- 
fine our discussion to a structure of the type shown in 
Fig. 4.1. Because the system possesses symmetric mass and 
stiffness properties, in addition to symmetrical boundary 
conditions, the solution of the eigenvalue problem consists 
of eigenvectors of two types, namely, symmetric and antisym- 
metric with respect to the center lines of the system. 
These center lines are chosen to coincide with the reference 
axes x and y, as shown in Fig. 4.1. It will prove conve- 
nient to work with a triple index notation (i,j,k) where the 
first two indices i and j represent nodal coordinates, as 
shown in Fig. 4.1 and the third index k represents the actu- 
ator type applied at the particular node. Moreover, k=l 
represents a force vector or translation in the z direction, 
perpendicular to the lattice plane, k=2 and k=3 represent 
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torque vectors or rotations in the x and y directions re- 
spectively. 

We recall that n & and n^ represent the numbers of beam 
members along the sides of the lattice and that n c is the 
number of cable members in each cable branch. With no loss 
of generality, we consider one element per member so that 
the order of the system in Eq.(2.57) is n=3 ( n a +l ) (n b +1 ) . 


The modal matrix can be written as 


U=(Ui u 2 . . .u . . .u ) 
~ ~c ~n 


(4.1) 


where 


u. 


f fa fb 

u 2. ! ~ 2 » “ 2 ’ 1 j 

f fa fb 

U «. [~ 2 * 2 


, 2 


n n, 

a b 


u * r 2 » ~ *> • 3 


n n, 

"i |- 2 s + 1 ' • T ’ 1 


(i.j.k) 


• m n, 

a b ~ 

U 2 . 2 ’2 ’ J 


2. 6 (1,2, . . . ,n} 

f n n 

i e l ' 

k 6 {1,2,3} 


~ - 1 


fb 

2 


— - 1 


(4.2) 

•£} 

■r} 
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in which the arguments of the entries in are ordered tri- 
ples of rationals (i,j,k), ordered in increasing lexico- 
graphical order of the ordered triples (j,i,k). 

Mote : (ix/ji/ki) preceedes (i 2 ,j 2 ,k 2 ) in lexicographical 

increasing order if and only if 
i i <i 2 

or, ii=i 2 and jj<j 2 

or, ii=i 2 and j t = j 2 and kt<k 2 

By virtue of the symmetry and antisymmetry of the modes, we 
introduce the following relations: 

Definition 

A mode is said to be symmetric if 
u(i,j,l)= u(-i,j,l) 

u(i,j,2)= u ( -i , j , 2 ) ( 4 . 3a, b, c ) 

u(i, j ,3)=-u(-i, j ,3) 

implying 


u(0, j , 3 )=0 (4.3d) 

A mode is said to be antisymmetric if 
u( i , j , 1 )=-u( -i , j , 1 ) 

u(i, j ,2)=-u(-i, j ,2) ( 4 . 4a, b, c ) 

u( i , j , 3 )= u( -i , j , 3 ) 


implying 
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u(0,j,l)= u( 0 , j , 2 )=0 ( 4 . 4d ) 

The matrix 

U T — [ u j u 2 ... .u ] T (4.5) 

C ** c 

contains the first c rows of the transposed modal matrix U. 
As far as the actuator placement matrix B is concerned, we 
let 

{ ( ii / j i / kj. ) / (iz* 3z* iz) * * • J ( ^ * ^ ) 

be a subset of the set of all arguments of the entries of 
any column in the modal matrix U. 


We recall the definition of the standard unit vector 


e n =(0,..., 0,1,0 Of (4.7) 

~ i 

which is a vector of n coordinates with a unit in its ith 
entry and zeros elsewhere, ie{l,...,nj. Using Eqs.(4.6) and 
(4.7), the actuators matrix can be written as 


B = [e n (ii,j! ,ki) 


e n (i ,j ,k )] 
> v c c c 


(4.8) 


B has full rank if no two ordered triples of E in Eq.(4.6) 
are identical. Multiplication of the matrices in Eqs.(4.6) 
and (4.8), results in the product 
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“l (il ,J1 ,ki) 

U l (i c’ J c’ k c ) 

Mil 

• 

• 

... Mi c ,j c ,k c ) 

u c (il,jl,k 1 ) 

— u (i ,j ,k ) 
c c J c c 


( 4 . 9 ) 


T 

We recognize that U C B is a minor of the transposed modal ma- 
trix U. We propose to use the shorthand notation 


T T 

u c B = u 


l(il.jl.kl) (i2»j2»k2> ••• (i c ’j c » k c ) '' 


( 4 . 10 ) 


Here, the first row indicates the selected rows of U T and 

the second row indicates the selected columns of U T in pro- 

T 

ducing the minor U B. 


The actuators placement over the domain can have either 
symmetric or asymmetric configuration. 

Definition 

We state that a set z of actuators exhibits a symmetric con- 
figuration if 

(i,j,k)eZ implies (-i,j,k)6H (4.11) 


Remark 
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Without constraints, any ordered triples of E of the form 
(0,j,k) do not affect the symmetry of the configuration of 
actuators . 

As illustrations, we selected six cases displayed in 
Fig. 4.2. The control scheme provides control of the first 
six modes by means of six thrust and torque actuators. 

Three cases represent asymmetric and three symmetric actua- 
tors configurations. Cases a, b, and c present a variety of 
asymmetric actuators placement configurations. It is not a 

difficult task to show by examining the determinant of the 

T 

matrix, that all the three cases yield regular matrices U^B. 
On the other hand, the last three cases d, e, and f present 
a variety of symmetric actuators placement configurations. 
Although cases d, e, and f are of the same kind, namely, 
symmetric, the nature of the resulting matrices U^B is not 
of the same kind... J2ases d and e yield zero determinant of 
the matrix U^B, proving U^B singular in both, while case f 
yields a non-zero determinant, proving U^B for this case re- 
gular. Cases e and f seem to be similar, although the first 
resulted in a singular matrix, while the second in a regular 
one. The question arises now as to how to distinguish beet- 
ween such two cases and how to avoid singularity. 

Notation 

Let E be a symmetric configuration of actuators 
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Notation: O force vector in the z direction (k=l) 

torque vector in the x direction (k=2) 

t torque vector in the y direction (k=3) 


Figure 4.2: Various Actuators Configurations 
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Figure 4.2: 


T 

Nature of U^B 


singular 


singular 


regular 
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- We denote by 2p the number of actuators in H of the form 
(i,j,k) such that i*0. 

Remark : There are exactly p actuators in H having i>0 

Remark : There are exactly c-2p actuators on the axis of 
symmetry i.e., the y axis. Note that c is the total number 
of actuators which is also the number of controlled modes. 
From this follows that there are in E exactly c-2p actuators 
(i,j,k) having i=0. 

- We denote by h the number of actuators in S of the forms 
(0,j,l) and (0,j,2), for some j. 

- We denote by v the number of actuators in 5 of the form 
(0, j , 3 ) , for some j. Clearly, c=2p+h+v 

- We denote by s the number of symmetric modes in the modal 
matrix U among the selected first c controlled modes. 

Proposition 

Let E be a symmetric configuration of actuators with dispo- 
sition 2p, h, and v. A sufficient condition for the mode 
. . . T 

participation matrix U C B to be singular is 

p+h<s (4.12) 

Explanation 

The proposition states that for the matrix U^B to be singu- 
lar, the number of actuator pairs p, in addition to the num- 
ber of actuators h with force or torque vectors, not paral- 
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lei to the axis of symmetry, must be less than the number s 
of symmetric modes controlled. 

Next, let us examine the cases of Fig. 4.2. The order 
of the system in Fig. 4.2 is n=45 . The solution to the ei- 
genvalue problem can show that among the first six eigenvec- 
tors, four have symmetric and two have antisymmetric modes 
shapes, so that s=4. A summary of the cases in Fig. 4.2, is 
given in Tab. 4.1. It is clear from the table that there is 
no conflict between the results obtained using the theorem 

and those obtained by inspecting directly the determinant of 
T 

the matrix U C B. 

Conjecture 

If z is not a symmetric configuration or if a is a symmetric 

T 

configuration with p+h£s, then U C B is regular. If this con- 
jecture is true, then the condition in the theorem is neces- 

T 

sary and sufficient for U C B to be singular provided that z 
has a symmetric configuration. Otherwise, if z is not a 
symmetric configuration set, the conjecture asserts that U^B 
is regular. Before proceeding with the proof of the theo- 
rem, we will mention the following important lemma 

Laplace Expansion of a Determinant 

The classical Laplace expansion for determinants has the 


form 
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TABLE 4.1 

Summary of Cases from Fig. 4.2 


Case 

p 

h 

V 

p+h 

s 

Result 

T 

u c B 

a 





4 

E is not , 

a symmet- 

H 





A 

ric configuration. 

D 






therefore 

p, h, and 







v are undefined 

c 





4 



d 

2 

1 

i 

3 

4 

p+h < s 

singular 

e 

2 

1 

i 

3 

4 

p+h < s 

singular 

f 

? 

2 

0 

4 

4 

p+h = s 

regular 
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r , w=l K K 

’il i 2 ... i' 


\ » i * ,* t ^ 

■1*1 lo ••• J_ 

1 L n-p 

I (-1) det A 

P 

*det A 

such chat p 

Jl 12 ••• j , 


••• K-? 


(4.13) 


where 


it 12 


Ui 3 2 


denotes the minor of the matrix A obtained by the selected 

rows i 1# i 2/ ...,i and selected columns j t , j 2 , . . . , j p . 

In the sum, there are ( n ) terms as the number of selection 

P 

of ordered p indices. The sum is extended over all p-tuples 
1£ j x < j 2 < . . . < j p^n, where j'i , j 2 / • • • 3^_ p denote the complemen- 
tary set of indices to j 1/ j 2 ,...,j p in 1,2, ... ,n arranged in 
natural order, and similarly i l , i 2 , . . . , i^ denote the com- 

plementary set of indices to i 1( i* i p in 1,2 n (see 

examples in Appendix A. ) 


Proof of Proposition 

Consider a set E with a symmetric configuration of actua- 
tors, with disposition 2p, h, and v. This represents a set 


E of the form 
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“ * L (ij » j 1 i^l ) t ••• 
(°»jp + i* k p+ i)» 
^’^p+h+l’ 3 ^ ’ * 


• ( VW’ 

•••* (0,:5 p+h’ k 

... (0,j c ,3)} 


( _i l »J1 . k l) . 

p+h )j 


( 'VW- 


(4.14) 


where are not zero, and k p+ x , . . . k p +h e { 1 , 2 } . 

T 

Let p+h<s hold true. We have now to show that U B is 

C 

T 

singular, i.e., that det(U B)=0. By virtue of the fact that 

B 

a change of rows (columns) in a determinant only affects the 
sign of the determinant, the vanishing of det(U^B) is invar- 
iant under such changes. Let us perform changes of rows in 
U^B, such that the first s rows will represent the s symme- 
tric modes and the rest c-s rows will represent the antisym- 
metric modes, renumbering the rows accordingly. 

T 

Let us now perform changes of columns in U C B such that 
the new order of the columns will be 

(ii - j i / k x ), (-ij , j j ,kj ),..., (i p , j ,k ) , (-i , j ,k ) 2p triples 

( 0 / j p+i / h p_j_i ),...,( 0 , j p+h / hp_|_h ) h triples 

(0,j p+h+1 ,3),...,(0,j ,3) v triples 

c 


Explicitly, 
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dec (Uj B) = 


block P 



u l(il»jl. k l) 

U 1 M-l ,Jl » k i) 

... Ui(i p ,j p ,k p ) 

Ul( ' i P ’ j P ’ k P ) 


u 2 (^1 » j 1 > k l ) 

u 2( _i l » J 1 » k l) 

U 2 (i p’ j p’ k p ) 

U2( ' i p’ j p’ k p ) 

a*det 

u s (il Jl , k l) 

u s (-ilJl> k l) 

••• wvv 

u (-1 ,j ,k ) 
s p’ J p* p 


u s+l (il *jl’ k l ) 

• 

• 

u g+1 (-il,jl, k l) 

••• u s+l (i p’W 

Vi ( 'WV 


• 

L« c (ii»Ji. k l) 


... u (i , j ,k ) 
c p p p 

V'VVV 


block H 

block 

V 

“i<°-WW •• 

• u l (0 -Vh’Vh ) 

1 Ul(0 ’ ;i p+h+l’ 3) 

• • • (0 » j ^ > 3) 

, U- (0 

1 

• 

i * 

,j p+l ,k p+l ) " 

• u 2(°»j p+h . k p+h ) 

1 U2(0, jp+h+l’ 3) 
1 

. . .a 2 (0,j c ,3) 

“s<° 

,j p+l’ k p+l ) " 

• U s <0 -Vh >l W 

1 U s (0,j p+h+l’ 3) 

. . .u (0,1,3) 

s c 


u j-i (0, j , i ,k , . ) u . , (0,i ,k . , ) 

s+1 J p+1* p+1 s+1 * J p+h* p+h 


* * * U C (0» j„j.u 


'p+h* p+h 


“s+l^ 0, Jp+h+l ,3) •' •“s+l^’^c’ 3 ' 


u c (0,:i p+h+l' 3) ••• u c ‘ 0>J c' 3) 


where o 6 {+1, -1} 


(4.15) 
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Note the blocks division P, H, and V in Eq.(4.15) where 
block P relates to the actuators located off the axis of 
symmetry, block H to the actuators type k=l and k=2 that are 
on this axis and block V relates to the actuators type k=3 
located along the axis of symmetry of the domain. 

T 

Noticing the fact that the vanishing of det(U c B) is in- 
variant under rows (columns) elementary transformations, we 
make use of the modes symmetry properties, as described by 
Eqs.(4.3) and (4.4) to perform such operations on the co- 
lumns of the determinant in Eq.(4.15). 


dec (Uj B) 


block ? 


= a*det 


••• u i ( vvV 

...u2Ci p .3 p .k p ) 


0 

0 


uidi.ji.kx) 0 

u 2di,ji,ki) 0 

UgCii.3i.ki) 0 

u s+l (il ’31 p l 2u s+l (i l »3 1 *kl ) • • - U s+1 ( W V P p 2u s+l ( V j p’V 

# 

« 

[u (i lt ji,ki) p 1 2u c (i 1 ,ji ,kx) • . * u c (ip’3 p .k p ) P p 2u c (l p ,:i p ,k p ) 


••• u s (i p ,j P’ k p ) 
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block 

_ 

H 



block V 


itn (o, 

i 

,j p+l’ 

VP 

• • • 

U 1 

(0 

’^p+h’kp+h^ 

|0 

1 

N 

0 

lu 2 (0 

i • *" 

• 

3 P+1 

k p4-l ) 

• • • 
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'^p+h’ k p+h^ 

1 0 • • • 
1 

0 

i* 

lu (0, 
1 3 

,j p+r 

•VP 

• • • 

u 

s 

(0 

’ ^ p+h ’ k p+h ^ 

1 

1 0 • • • 

1 

0 

10 

1 • 




0 



[ u s+l (0,j p+h+l’ 3) 

u s+l (0,:i c ,3) 

I* 

10 



• • • 

0 



! 

|u c ( 0 ‘ j p+ h +r 3) ••• 

u c (0J c ,3) . 


where cr € {+1, -1} (4.16) 


The final form as obtained in Eq.(4.16) is a result of the 
following operations: 

Consider column 2 2, and column 21-1 for any 2 e (1,2, . . . ,pj . 

If k is 1 or 2, then in order to achieve zeros in the first 

2 

s rows in column 2i of block P, we subtract column 22,-1 from 
column 22,. If k^ is 3, then in order to achieve the same 
goal, we add column 22.-1 to column 2 . The columns 22,-1 for 
2 . 6(1,2, . . .p} remain as they are. , The parameter that ap- 

pears in block P, represents a sign function 
p (k) = (-l)^ 3-k ) such that p (l)=p (2)=-l and p (3)=+l. 

5, A* Ac X 

Applying Laplace expansion Eq.(4.13) to the determinant of 


Eq. (4.16) 
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det(U^B) 

= a*det(U^B) 

= 

o:(-i) 1+2+ 

. . +S+j 1 +J 2 + 

...+3 S 



/ 

I 




T 

1 

2 . . . s 

m 

s+1 s+2 ... c 


• dec U^B 



• det U^B 

• » • f . i 



3 1 

32 ••• 3 S< 




where the sum is as indicated in Eg. (4.13). In this expan- 
sion, U^B denotes the rearranged matrix as in Eq.(4.16), the 
rows (columns) numbering is also with respect to the new ar- 
rangement, and the choice of first fixed s rows was made. 

Claim 


det U^B 


'l 2 
Ji iz 



(4.18) 


for any j 1 ,j 2 ,... / j s such that 1£ j x < j 2 < . . . < j 

Proof of Claim 
The matrix 

1 2 ...s' 

U t B 

C 

U 

has at most p+h nonzero columns. By assumption, p+h<s. 


Therefore, this matrix of order s has at least s-p-h>0 null 
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columns. Because s y p, and h are integers, s-p-h is an in- 
teger also, and a positive one. Hence, at least one column 
of the above mentioned matrix is zero, so that the determi- 
nant of the matrix is zero. This concludes the proof of the 
claim. 

In view of this claim, we have 

det(U^B)=0 (4.19) 

because its expansion in Eq.(4.17) is really a sum of zeros. 

T 

This proves that U^B is singular, which concludes the proof 
of proposition. 



Chapter V 

MINIMUM FUEL IN HIGH ORDER SYSTEMS 


The minimum fuel problem became a very important one in 
recent space vehicles, especially in those that are launched 
for lengthy missions. The amount of fuel or energy alloted 
for the control system is limited to such a degree that it 
becomes necessary to treat fuel economy as the predominant 
factor in many systems. In such cases, it is natural to de- 
sign the control so that it consumes a minimum amount of 
fuel . 


There are many problems where the control variables 
F(t) are directly proportional to the rate of fuel flow, 
like the flow from a gas jet for attitude control or shape 
control of a large flexible spacecraft. Almost invariably 
in these problems, the total fuel available is limited. 
Therefore, it is desirable to accomplish each control cor- 
rection with a minimum amount of fuel. The total fuel con- 
sumed during control action is measured by the time integral 
of the absolute value of the control variable. Due to tech- 
nical restrictions on the actuators, there are limits on the 
magnitude of the control forces. As an example, the thrust 
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produced by a gas jet actuator is limited in magnitude due 
to the saturation of the power elements. Systems character- 
ized by a restriction on the control and subjected to the 
minimization of a particular functional, such as the con- 
sumed fuel, the required energy, etc., are optimal. We re- 
fer to a control system minimizing the amount of consumed 
fuel, as a "minimum-fuel system". 

About two decades have passed since the development of 
powerful variational techniques in the theory of control, 
such as the minimum principle of Pontryagin. This principle 
can be used to determine necessary and sufficient conditions 
for the control function to be optimal [16,17], The neces- 
sary condition usually specifies the nature of the control 
and the general structure of the control system. Unfortu- 
nately, the derivation of the optimal fedback control law, 
i.e., the explicit dependence of the control on the instan- 
taneous values of the state variables, is a very complex 
task and for most practical cases it remains an unsolved 
problem [18,19]. In these cases, a trial and error process 
is unavoidable for the state determination. One may consid- 
er such a process as reasonable for a low-order system 
(fourth order at most), but certainly not for a high-order 
system, where the computational difficulty is insurmontable . 
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However, it is unexpected that a reasonably accurate model 
of a flexible structure can be represented by a low-order 
system. 

Formulation of the control law in terms of the costate 
variables can be accomplished regardless of the system ord- 
er. Yet, determination of the initial conditions of the cos- 
tate variables for a high-order system is difficult, if not 
impossible. The control task is made considerably simpler 
by using the IMSC method, where the complexity inherent in a 
high-order system is reduced. Using this method, the ini- 
tial constate variables must be determined for a set of in- 
dependent second-order ordinary differential equations, so 
that treatment of the coupled high-order system is avoided. 
According to this approach, to each of the decoupled modes, 
independent control policy can be specified so as to permit 
independent control of each individ a-L- mnde — With this meth- 
od, the cumbersome minimum fuel problem in the high-order 
system is reduced to n fuel minimization problems of the de- 
coupled set of n second-order modal space systems. 

After the control design task has been carried out in 
the modal space, and modal control functions have been 
found, these generalized control forces are synthesized to 
obtain the actual actuators forces. 
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5.1 PROBLEM FORMULATION 

In formulating the minimum fuel problem, we begin with 
the discretized time- invariant structure equations of motion 
of order n 

M q(t) + K q(t) = Q(t) (5 1} 

where Q(t) is the discrete forcing vector, generated using c 
actuators, or as many as the number of controlled modes, 

Q(t) = B*F(t) (5.2) 

where F(t) is the c-dimensional vector of thrusters and 
tourquers and B is a full rank nxc actuators placement 
transformation matrix. 


The problem of minimum-fuel control can be stated as 
follows: Assume that the control F^(t) is proportional to 
the rate of flow of fuel. Then the total fuel consumed J(F) 
during a time interval (0,t f ) that is also a measure of con- 
trol effort, can be expressed in terms of the time integral 
of the sum of the absolute control values as follows: 


ft f 

( \ 
C 

J(F) * 

J 0 

l ajF (t)| 

i=l 1 1 

V 4 


(5.3) 
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where are nonegative weighting factors and tf is some 
unspecified final time. In most practical systems, the mag- 
nitude of each of the control functions F ± ( t ) ( i=l , 2 , . . . , c ) 
is limited due to the saturation of the power elements. 
Therefore, it is assumed that the controls have to satisfy 
the constraints 


|F (t)| < M i - 1,2, ...,c (5.4) 
where M are positive constants. 

It is assumed that at t 0 =O the system (5.1) is in the 
initial state q(t 0 ), q(t fl ). Then, the control task amounts 
to driving the system from its initial state to a final spe- 
cified state q(t^), q(t^) in finite time. 

The essence of the control problem is to determine a control 
vector F(t) subjected to the constraints (5.4), such that 
the final state is reached with a minimum amount of fuel 

J(?). 

As stated already, the problem can be decomposed into a set 
of n second-order decoupled optimal control problems, which 
reduces the amount of complexity arising in high-order sys- 
tem. In the formulation, we employ theoretical results de- 
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veloped in Refs. 16,17 to obtain the necessary condition for 
the optimal minimum- fuel control. 


The modal linear time-invariant plant is described by 
the matrix differential equation 


x (c) = A x (t) + b N (t) 
~r r -r - r 


where 


r = 1,2, ...,c (5.5) 


x 

-r 




9 


A 


r 


0 



1 

0 


b = [0 1] T 


(5 . 6a, b, c) 


x lr (t) 


5 r (t) 


2r 


e r (t) 


(5 . 6d) 


5 r (t) and £ r (t) are modal displacement and modal velocity 
respectively, x r (t) is a 2x1 vector, A r is a 2x2 matrix and 
N r (t) is the scalar modal control function. 


To apply the IMSC method, one must seek the minimiza- 
tion of modal cost functions rather than the minimization of 
the global cost given in Eq.(5.3). To derive independent 
modal cost functions, a decomposition of the global cost 
into modal costs, dependent on the associated modal control 
forces and modal coordinates only, is essential. 
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It turns out that a complete decomposition of the global 
cost for fuel minimization problem, as it appears in 
Eq . (5.3), into equivalent independent modal costs is not 
feasible, so that we resort to another approach. 


Using Eqs. (2.72a) and (5.2), the modal control vector 
is obtained as 


N c (t) = C F( t) (5.7) 

where 

C = U*B (5.8) 


is the system cxc mode participation matrix. We employ 
norm-one measurement for the configuration space actuator 
force vector F(t). By definition, norm-one is the sum of 
the absolute values of the vector entries which resembles 
exactly the global cost presented by Eq.(5.3), or 


jF(t)j = l If (t) 


(5.9) 


i=l 


This norm is defined for the n dimensional Euclidean vector 
space R n and has all the properties of a regular vector 
norm, as shown in Ref. 29. 

The related matrix norm is defined as 


II C || = max y c . . 


j i=l 




(column sum) 


(5.10) 
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Using these norm definitions we can write the norm-one 
relation, applied to the control vector transformation in 
Eq . ( 5 . 7 ) , as 

|F(t)| < IlC- 1 !! |N c (t)| (5.11) 

With no loss of generality, we assign the values a, =1 for 
all i=l,2,...,c as weighting factors in the cost Eq(5.3). 
Note that the derivation can be easily modified if those as- 
sumptions are not made. We propose to minimize the modal 
fuel expenditure by minimizing the norm-one of the modal 
control vector. This vector consists of independent control 
entries, so that minimization of its norm-one can be carried 
out by minimizing each entry independently. Integrating 
Eq.(5.11) on (O/t^) will result in a cost inequality that 
provides an upper bound to the global cost of Eq.(5.3) 



its modal counterpart. 

It was stated earlier that the magnitude of each actua- 
tor force is bounded due to limitations on the forces char- 
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acteri sties, as can be seen in Eg. (5.4). These bounds must 
be transformed to modal forces bounds before the modal mini- 
mum-fuel problem can be formulated. We use the norm of 
Eg. (5.7) to determine modal control bounds. Eg. (5.4) in 
conjunction with the norm-one definition of Eg. (5.9) provide 
the bound to the norm of the actuators force vector, or 


FCOj - I |F (t)| < i H 
i=l i-1 


(5.13) 


Applying norm-one measurement to Eg. (5.7), we obtain 


N (t)| < IICII | F(t) | < IICII l M 
~ c ' i-1 


(5.14) 


and, because |N ( t ) \ x =1 !n . ( t )| , it turns out that 

i=l 


l |N.(t)| < IICII V M 
i-1 1 i=l 


(5.15) 


which is a total bound. 

In the IMSC method, the modal controls act independently of 
each other, so that each control must have its own bound. 
To achieve this requirement, we introduce weighting coeffi- 
cients a such that 

i 


C 

IN. I < a. IICII l M. 

J-l J 


(5.16) 


where 


a. > 0 

i 
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One way of generating these coefficients is to select them 
inversely proportional in some manner to the eigenvalues. 
In a case where n r rigid body modes are present, as well as 
c-n r elastic modes, one may choose 


a. 

i 



1,2,... »n^ 


(5.18) 


for the rigid body modes, and 


a. 

x 


r 


(c-n r ) 

c(c-n r -l) 


X. 

l 


c 



i = n r +l,n r +2,...,c (5.19) 


for the elastic modes. 

Clearly, any other set of weighting coefficients can be used 
as long as Eq.(5,17) is satisfied. 


The essence of the IMSC method is to design the control 
laws in the modal space for each of the controlled modes in- 
dependently. Then, realization of the actual force vector 
is made by modal synthesis. It can be shown [17,18] that 
for modal minimum-fuel problem the controls appear as a non- 
linear "bang-off-bang" type, that use the maximum values of 
the modal control forces. This means that the norm of the 
modal control vector will reach its upper bound in Eq.(5.15) 
very often. In such case, through substitution of Eq.(5.14) 
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into Eq.(5.11), we see that the norm of actuators forces 
would exceed its bounds, or 


F(t) | < II C -1 !i HCII l M. 

i=l 


(5.20) 


— 1 . 


where ||C - ||||C|| >1 is the condition number of the mode partici- 
pation matrix C, which is unacceptable. To circumvent this 
problem, it is necessary to modify Eq.(5.13) by dividing 
each control bound by the condition number of the matrix C, 
so that 


|F(t)| < 

lie ii ti cu 



(5.21) 


Eq.(5.21) provides the effective upper bounds on actuators 
forces, where its modal counterpart will be obtained as the 
modified Eq.(5.14) 


|N C (0| 




1 


lie 




(5.22) 


Clearly, the condition number plays the role of a "figure of- 
merit" in the control effectiveness. 
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Now, we are ready to formulate the modal minimum-fuel 
problem for the linear time invariant modal plant described 
by the modal differential equations (5.5). Because of limi- 
tations on actuators power, the norm of the modal control 
force vector is bounded as given in Eq.(5.22). Using the 
weighted decomposition scheme given by Eqs . (5.16-19), the 
modal forces bounds can be given for each mode as 


c 

|N (t)| < a — l M. = N r = 1,2,. ..,c (5.23) 

r r IIC" II i=l 1 r 

The task of the modal controls is to drive the initial modal 
states x 0r to final desired states (the origin for exam- 

ple), by appropriate optimal control signals N r (t), minimiz- 
ing the performance indices J r (N r ). In our case, the modal 
fuel functionals are given by 


J (N ) = 
r r 


N r (t)|dt 


r = 1,2, ...,c (5.24) 


The minimum principle of Pontryagin [16] can be used to der- 
ive the necessary condition for the optimal modal minimum- 
fuel control. Using Eqs. (5.5) and (5.24), we introduce the 
modal Hamiltonian functional 
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H r (x r ,p r ,N r ,t) = |N r (t)| + p*(t) (A r x r (t) + b N t (t)) (5.25) 


r * 1,2, ... ,c 


where p r ( t) = [p lr ( t) p 2r (t)] T is the 2x1 vector of Lagrange 
multipliers known as the modal costate vector. 

The necessary conditions for minimization are given by 


x r (t) = 

T 9p, 


c e (0,c f } , r = 1,2,. ...c (5.26) 


P r (t) 


c 6 (0,t f > , r = 1,2, ...,c (5.27) 


c 6 , r = 1,2, ...,c (5.28) 


subject to the boundary conditions 


x (0) = x 
-r -Or 


(5.29a) 


x (t c ) = x. 
-r f -fr 


(5.29b) 


H r ( -r <t f ) ’?r (C f ) ’ N r (t f ) ’ t f ) " 0 


r = 1,2,. ..,c (5.29c) 


The object is to determine the 2c constants of integration 


and t 
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Because the cost functional is not a function of p r (t), 
it follows that Eqs.(5.26) are identical to Ecs.(5.5), the 
modal plant equations. Eqs.(5.27) lead to the adjoint-cos- 
tate system of differential equations 


?r (t) = ' A r ?r (E) 


r = 1,2, ...,c (5.30) 


The control N r (t) minimizing the Hamiltonian H r absolutely 
is a function of x r , p r and t. This functional dependence 
of N r is a necessary condition on the control function 
N r (t), which is the one that minimizes the performance index 
J r^ N P' can be stated that the optimal control N r (t) 
minimizing the performance index J r (N r ), Eq.(5.24), must by 
necessity minimize the Hamiltonian H Eq.(5.25), absolute- 
ly- 


The formal minimization technique, Eqs.(5.28), leads to 
an unrealistic result. However, 



r = 1,2,. ..,c (5.31) 

does provide an absolute minimum to the Hamiltonian func- 


tion. 
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The shaded area in Fig. 5.1 provides the value of 
| N ( t ) | +p^(t) bN r ( t ) versus p r (t)b depending on the value of 
N r ( t ) . From the plot, it is clear that H ^ takes its small- 
est value at each instant t if the term |N r ( t ) | +p^ ( t )bN r ( t ) 
maintains its largest negative value at all times. This can 
be achieved by selecting N f (t) along the solid lines provid- 
ing the lower bounds of the shaded area. Hence, we consider 
the switching functions 

N r (t) - - N r sgn (p*(c) • b} Ig^(t) b| > 1 (5.32a) 

V C > = 0 l?r (t) < 0 (5.32b) 
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Figure 5.1: N-Dependent Part of Modal Hamiltonian Function 
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An illustration of Eqs.(5.33) for the rth mode fuel-optimal 
control is given in Fig. 5.2. Equations (5.32) provide the 
modal control law for the modal minimum-fuel operation, whe- 
re N r (t), (r=l,2, . . . ,c) are piecewise constants functions 

A A 

with values of +N r ,0, or -N r , as given in Eqs . ( 5 . 32a, b ) . If 
however, there are time intervals (t^tjJr of finite dura- 
tions, during which one or more modal necessary conditions, 
given by Eqs. (5.31), provide no information about the behav- 
ior of the controls N r (t), then are singular inter - 

vals and the appropriate N r (t) are singular controls in 

these particular modes. Such a situation occurs when 
T 

I p r b |=1. We observe that Eqs. (5 . 32a, b) specify the magni- 
tude and the polarity of the fuel optimal control N r (t) in 

T T 

terms of p r (t) uniquely, provided Ip^blxl. When |p r b|=l, as 
in Eqs . ( 5 . 32c , d) , the fuel-optimal control polarity is spe- 
cified, but not the magnitude. These considerations lead to 
the distinction between normal modal fuel-optimal problem 
and singular modal fuel-optimal problem. 

We use results from Ref. 16 to show that singular in- 
tervals can exist only in the controls of rigid body modes 

but never in the controls of elastic modes. The condition 

T 

is based upon the matrix E A where A is the modal plant 

r r r 

matrix and E^_ is the 2x2 modal controllability matrix, given 
as 




Figure 5.2: Modal Fuel-Optimal Control Function 


i cr 
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E = [b i A b] ( 5 . 34) 

r i r J 

It is shown, that a necessary condition for a singular in- 

T T 

terval to exist requires that the matrix [b A r b] A be sin- 
gular. This occurs when either the modal system, Eqs.(5.5) 
and (5.6), is not completely controllable, in which case 
(b A^_b] is singular, or when the modal plant matrix A f is 
singular. The modal controllability matrix has the non- 
singular form 


E = 
r 


0 

1 


1 

0 


(5.35) 


for both rigid or elastic modes, provided all modes are con- 
trollable. In contrast to the controllability matrix, the 
modal plant matrix A r , given by Eq.(5.6b), becomes singular 
in the case of a rigid body mode for which u> r =0. 

Hence, a necessary condition for a singular interval to ex- 
ist is for the system given by Eqs.(5.5) and (5. 6) not to be 
completely controllable, or for the plant matrix A r to be 
singular. 

In view of the above discussion, it is clear that each 
of the two mode-types, namely, rigid body modes and elastic 
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modes, possesses different properties. Therefore, it is ap- 
propriate to discuss the design of modal fuel-optimal con- 
trol policies for these two unlike modes separately. 


5.2 FUEL-OPTIMAL CONTROL FOR A RIGID BODY MODE 


The control of a rigid body mode (RBM) in the modal 
space, resembles a single-axis attitude control of a unit- 
inertia space vehicle. Rigid body modes from the set of de- 
coupled differential equations for the number of controlled 
modes described by Eqs. (2.77a), have the form 


5 ? (t) - N f (t) r » 1,2, ...,n r (5.36) 

where n is the number of rigid body modes. Using the same 
state variable definition as in Eq.(5.6d), the modal state 
vector x r (t) satisfies the matrix differential equation 


r. , i 

r- 


1 


) 

x ( t) 


0 1 


x (t) 

It* 


0 

I L 

► = 


< 


■ + - 

> 

x (e) 


0 0 


x„ (t) 


1 

2r J 

1 



2r J 


W i 


N r (t) r = 1,2, ...,n r (5.37) 


It is desired to drive the modal plant from an arbitrary in- 
itial state given by 


*lr (0) = ? r (0) = ? 0r 


x 2r (°) 


£ r (0) 


*0r 


1 , 2 ,. 


(5.38) 
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to a desired final state, say the origin, using a minimum 
amount of modal fuel, or 


J (N ) 
r r 


rt f 

| N (t) | dt 

J 0 


r = 1,2,. ..,n r (5.39) 


subject to the constraints 
|N r (t)| < N r 


r * 1,2, . . . ,n r (5.40) 


where N is the saturation value for the modal control 
r 

force, as defined by Eq.(5.23), and t f is the final time, 
which is free. 

To this end, introduce the modal Hamiltonian 


tf r (x r ,p r ,N r ,t) = |N r (t)| + p lr (t) x 2r (t) + p 2r (t) N r (t) 

r => 1,2,,. 


where, the costate variables p. (i=l,2) are the 

lr 

of the adjoint system differential equations 


(5.41) 


solutions 


r. i 



f 

p lr (t) 


o 

o 


P, (c) 

! 1 r 

II 



lr 1 

p 2 r (t) . 


-1 0 


P 2r (t) 


r = 1,2, . . . ,n r (5.42) 


Note that even if the controls N r (t) are given as a function 

of time, Eqs . ( 5 . 42 ) will not uniquely define a vector p (t), 

~ r 
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because both the initial and terminal (final) conditions 
have been imposed on x r< and no p r (0) boundary condition has 
been specified. 

The minimization of the Hamiltonian Hr with respect to 
N r (t), leads to the control law 

N^Ct) = - N r dez (p 2r (t)} r = 1,2, ...,n r (5.43) 

The optimization problem Eqs.(5.37) and (5.39) has thus been 
transformed into the problem of finding the solutions to the 
set of simultaneous Eqs.(5.37) and (5.42), with the control 
N (t) given by Eqs.(5.43), with the boundary conditions 

~r (0) = ? 0 r ’ W = 9 (5.44) 

specified, but with the final time tf unspecified. The 
knowledge that the optimum N f (t) has the form of Eqs.(5.43) 
simplifies the optimization problem immeasurably. An impor- 
tant comment should be made at this stage: The Pontryagin 

theorem provides a necessary condition for the optimal con- 
trol but not a sufficient one. In other words, if the theo- 
rem provides as a solution only one control function, this 
function is an optimum. However, more than one control 
function may satisfy the theorem, in which case it is neces- 
sary to search further for the optimum. 

The solution to the costate equations, Eqs.(5.42) becomes 
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P lr (c) - V (5.45a) 

P 2r (0 » n 2r - n lr C (5.45b) 

where 


V 


P, r <0) 


2r 


p._(0) 


r = 1,2, . . . ,n r (5.45c) 


are unspecified constants. 

Even though p 2r (t) is not known as a function of time, the 
mere fact that p 2r (t) is linear in time requires that N r (t) 
from Eqs.(5.43), must proceed in time as a nonrepeating se- 

A a 

quence of values of the form ( + N , 0,+ N } or any part of it. 
In view of this realization, if ii lr *0 the only nine modal 
control sequences that can be candidates for the fuel-opti- 
mal control are 


(0), (+N r ), {-N r >, {+N r ,0}, {-N^.O}, (0,+tO, {0,-N^} , 
{+N r ,0,-N r }, {-N r ,0,+N r } r = 1,2, 


,n r (5.46) 


Figure 5.3 provides a description of these control sequences 
versus their costate counterparts. A control sequence like 

A A A 

{-N r ,0,+N r } means that first N (t)=-N r , then N r (t)=0 and fi- 
nally N r (t)=+N r are applied in this order. Therefore, in 
order to transfer the state to the origin from any point in 
the state plane, it will be necessary to use one of the con- 
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trol sequences given in Eqs.(5.46). Clearly, any other con- 
trol sequence will violate the necessary condition for mini- 
mum- fuel control. 


Solution to the state equations is obtained simply 
through a direct integration of Eqs.(5.37) with the initial 
state given by Eqs.(5.38) for the optimal control given by 
Eqs.(5.43). The solutions are 


V (t) =l Cr + 5 or C 


*2r (t) ’ ’Or 


r » 1,2,. ..,n r (5.47a) 


r - 1,2 n r (5.47b) 


for N r (t)=0 and 


= 5 0r + «0r C + I V* 


*2r (t> * ’or + 4 r c 


for N r (t)=±N r =A r . 


r - 1,2,. . . ,n (5.48a) 
r = l,2,...,n (5.48b) 


State plane trajectories are then obtained by eliminating 
the time t from Eqs.(5.47) and (5.48) 

x, r = ^ 0r “ const , r = 1,2, ...,n r (^- -9) 


x 


lr 


- 5 


Or 


1 

2A 


5 0r 


1 

2A 


2 r 


r - 1,2, ...,n r ( 5 . 5 Q) 



Figure 5.3: Candidates for 



Fuel-Optimal Control for RBM 
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Equations (5.49) and (5.50), as well as Fig. 5.4a,b display 
these trajectories in the cases for which N r (t)=0 and 
N r (t)=±N r , respectively. Clearly, from Fig. 5.4, the only 
state that may be driven to the origin, using the control 
N r (t)=0 is the state (0,0) itself. Moreover, for any other 
state, it is necessary that either N r (t)=+N r or N r (t)=-N r 
will be the final control portion, approaching the origin 
along the trajectories labeled as Y + and T* in Fig. 5.4b re- 
spectively. 

From Eq.(5.50) it is clear that 


Y r ' ((x lr-' t 2r ) : *lr : X 2r * 0) 

r 


y = { (x, ,x„ ) : x, = - 
r lr’ 2r Ir 


1 ,2 


2N 


x~ ; x > 0 } 
2r 2r 


1.2.. .. ,n^ (5.51) 

(5.52) 

1.2. . . . ,n r 


The union of these two curves constructs the final switch 
curve for the modal minimum-fuel operation 


[ <x ir’ x 2r > 


ir 




r =* 1,2, . . . ,n r 


r 

Based on these observations, clearly the fuel optimal 
tion depends on the modal initial state (Sor/5or)- 
we observe four regions in the modal stare plane 


(5.53) 

solu- 

Hence, 




' r 

Figure 5.4: State Plane Trajectories for RBM 

(a) N r (t) = 0 A 

(b) N r (t) = ± N r 
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*1 - <<* 1 r-* 2 r ) : X 2r * 0 ! x lr > ' X 2r } 

ZN 

r 


(5.54a) 


R? = { (x ,x ) : x„ >0 ; x, < - -4 — xj; } 

*- lr’ 2r 2r ’ lr ~ 2r 

r 


(5.54b) 


■>3 - <<*, r ,* 2r ) ■ x 2r <» i * lr < 4-*t r > 


(5.54c) 


*“ ' l(x ir' x 2r ) : x 2r ' 0 ! X !r S X 2r> 


( 5 . 54d) 


These four regions are displayed in Fig. 5.5, as well as 
four cases of modal initial states, one in each of these re- 
gions . 

It is shown in Appen. B that if ( £o r , £ 0r ) e R i UR 3 ' then mini- 

■ 

mum-fuel control does not exist. If (Co r /5or ) 6R 2 UR 4 , then 
the control sequence {0,-N r } is a minimum-fuel control for 
(Sor,««r)eR 2 , and {0,+N r j is a minimum-fuel control for 
Uor,Sor)eR 4 . Hence, depending on the initial state, there 
may or may not exist a fuel-optimal solution. 


Fuel optimal systems with unspecified final time have 
the unfortunate property that the control resulting for the 
minimum expenditure of fuel, usually requires infinite final 
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time (See Appen. B) . We could accept limitations of this 
type if the space under consideration is the configuration 
space, by limiting the initial state of the problem with 
unspecified final time just to a state subspace R 2 UR 4 . In 
general, it is very unlikely to expect the modal initial 
conditions obtained by the linear transformation Eq.(2.66) 
to have a specified distribution in the modal space and 
still deal with a practical problem. Solution to the mini- 
mum-fuel problem must be capable of handling any initial 
conditions in the modal state plane. We can clearly con- 
clude from Eq.(B1.7) that, to assure the existence of a so- 
lution in region R t UR 3 , the response time, i.e., the final 
time t f , must be either fixed a priori, or bounded from 
above by a fixed time T f such that t f £T f . Thus we have to 
recognize the fact that bounding the response time will pe- 
nalize fuel expenditure. 

From the minimum-fuel solution, we can see that mode r 

arrives at the origin at a final time t f , depending on the 

initial state, such that t <t <», where t, 

f min min fr f mm mm 

is the lower bound, i.e., the shortest time period in the 
minimum-time solution to reach the origin among all cont- 
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rolled modes, and tf is the final time of the rth mode (See 
Appen. C for the minimum time solution for a second-order 
system). It is understood that the "modal mission" is in- 
complete as long as the last mode in the control scheme has 
not reached the specified target. Thus, the control task is 
to drive all the modes in the control scheme to the target 
(origin) in a bounded period of time. For the same reason, 
it is clear that no advantage is gained in driving one mode 
in a shorter time period then the others, consuming uneces- 
sary fuel. Therefore, we base our discussion on the premise 
that all the modes in the control scheme are designed to ar - 
rive at the origin simultaneously . 

Definition 

We shall say that T f is a minimum upper bound for the con- 
trol operation if 

C £ mil-max « T f < “ (5-55) 

where t' c is the largest minimum-time solution to 

f mm max 

reach the origin among all controlled modes. Certainly, T f 
may be selected larger then the minimum upper bound to meet 
some fuel constrains. 

Figure 5.6 displays a typical optimal path in the state 
plane for fixed final time with the control sequence 
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{-N r ,0 / +N r j. To specify the control uniquely, two switch 

times t and t must be determined. This is done in Ap- 
pen. C, where the switch times are obtained in the form 



(5.56a) 



(5.56b) 


From Eqs.(5.56) it is clear that both switch times depend on 
the modal initial conditions (Zor'Zor) as well as on the 
specified final time T f . 


Before describing the optimal control for the minimum- 
fuel problem with specified final time, we recall one of the 
results obtained in Appen. B. It has been shown in case 
B.2, for the initial modal state (Si* , 5 or )6R.4 / that the fi- 
nal time for the minimum-fuel problem, given by Eq.(B2.2), 
is 

(5.57) 





Figure 5.6: Modal State Plane Trajectories for RBM with Fixed Final 

Time 
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If the specified final time is Tf>tf, where tf is given by 
Eg. (5.57), clearly the state will reach the origin in less 
than the specified final time Tf with the minimum-fuel opti- 
mal control sequence [0,+N r j. In this case, an exception is 
made by permitting the state to arrive at the origin not at 
the specified final time, because any adjustment to the fi- 
nal time will either result in the same fuel consumption or 
it will increase the fuel consumption, in contradiction to 
the basic approach for fuel minimization. 

Related to the four state space regions R : to R 4 , the 
minimum-fuel optimal control law has the following expres- 
sions : 
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N f (t) =■ - N r for t 0 < t < O(e 0r ,? 0r ) € R X for all T f 

N r (t) =0 for t z < t < tWc 0r ,l 0r ) 6 R 4 for T f < t° 

N r (t) = + N f for c w < t < T f 

and we will select 

N r (t) = 0 for t 0 < t < c ll (5 0r *£ Qr ) 6 ^ for T f > c f 


N (t) = + N r for ti < t < t° 
N (t) =0 for t° < t < T c 


c i “ fc o + — 5. 


t 0 = tn _ _±_ | _ _or 

c f c ° ~ ^or 4 

2N r 0r «i0r 


’z 2 f 


T f + t 0 + 


T- - to - t 2 

X N 


Nr 2l 0r 


'w 2 1 f 


7 i T f + c 0 


I l 2 

„ ^ 0 r 

T- “ t(j ~ » 

N 


4£ 2E 2 1 & 

Sr Sr 





115 


N (t) » + N for tg < t < t ] 

r ^ Z 


N r (t) * 0 


for t < t < t 
z w 


(5 0r^0r ) 6 R3 for a11 T f 

l’ (5 0r^0r ) 6 R 2 for T f < c f 


N (t) = - N for t < t < T, 
r r w f 

and we will select 

N r (t) = 0 for t 0 < t < tj] 

.0 
•f 


N r (t) = - for tj < t < t° 


N r (t) = 0 
where 


for t° < t < T f 


U 0r'^0r ) 6 R 2 for T f > C f 


ti = tn - 


2N 


’Or 


Or 

^ 0 r 


l (5.58b) 


t0 f 


t 0 + 4- 5 - > 

2N r 0r 5 0r 


1 

2 


T f + c 0 ~ 


Or 


N 


T„ - 


to + 


’Or 


N 


4 C 


Or 


2| 2 


Or 


N 




r J 


C w=f < 


T f + t 0 - + 

N 


T„ - 


to + ~ 
N 


or 


4 € 


Or 


N 


2i 2 

_0r| 

N 2 

r 


It is shown in Appen. C that the locus of the first optimal 
switch time t z is a section of a second-order parabola in 
the forth quadrant, having the expression 



~ 

f 

• 

- 

Zi - - 

2N 

r 

T f ^ 0r 

T f c 0 - - — 

N 

r J 

Z 2 


(5.59) 
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For various values of the modal initial velocity, a family 
of parabolas will be necessary to describe the state at the 
first optimal switch time. Clearly, the second optimal 
switch time occurs when the state trajectory intercepts 
the I curve . 


Definition 

An initial modal state is said to be in reachable zone in 
the state plane, if the origin can be reached using the ap- 
propriate optimal control sequence in the prescribed final 
time bound. 


It is shown in Appen. C that the reachable zone is bounded 
by two sectors of a second-order parabola 


! • - 4 - [z| + 2(T f - t 0 )M Z 2 - (T f - t 0 ) 2 

4N L r J 


(5.60) 


and its reflected image about the origin. To guarantee that 
all rigid body modes are inside the reachable zone, the fi- 
nal time can be determined by the modified Eq.(5.60), or 


> to + max 


l N r L N N j J 


1^1,2, . . . ,n 


(5.61) 


for a modal initial state (So r ,So r ) eR i UR « • Similarly, for 
(Cor /5or )eR 2 UR 3 -, we obtain 


Tf ^ t 0 + max 


c— 1,2, . . . ,n 


{.^ + r 2 

\ N r L In N J ' 

2 n ^ r ' r 


(5.62) 
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Clearly, Eqs.(5.61) and (5.62) resemble the minimum- time 
solution, Eqs.(C1.5) and (Cl. 6), respectively. 

A dimensionless plot of Eqs.(5.51), (5.57), (5.59) and 
(.5.60) is shown in Appen. C, Fig. C.2. The plot provides 
also a description of reachable zones for the various con- 
trol regions, and switching curves for che optimal control 
law, Eqs. (5.58) . 


5.2.1 Fuel Consumed - RBM 


The fuel consumption depends on the initial condition 
point E in the modal state plane. Refering to Fig. C.2, for 
a point in region R' (R + ), the modal fuel depends also on 
the first and second switch times tj. and t w as shown in 
Fig. 5.6 and is given by 


J (N ) 
r r 


N dt + 
r 


r c f t 


N dt = N (t - t 0 + t, - t ) 
r r z 0 f w 


(5.63) 


For a point in region R° , the modal fuel is given by 


J 

r 


(N r ) 



N r (tO 


Cl) 


(5.64) 


Expressions for the times t fl , , t f , t t and t® are pro- 

vided in Eqs. (5.58). The total modal fuel for the RBM's is 
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obtained by summing all modal fuels in Eqs.(5.63) and 
(5.64), depending on the different modal initial conditions, 


or 

J(Nn r ) » l VV 

r=l 


(5.65) 


5 . 3 FUEL-OPTIMAL CONTROL FOR AN ELASTIC MODE 

Formulation of the control problem for an Elastic Mode 
(EM) has a similar pattern, to the one of a rigid body mode. 
The control of an elastic mode in the modal space is tanta- 
mount to controlling a set of linear oscillators system. 

The derivation of the control law in terms of the costate 
variables can be accomplished very easily. However, the 
task of determining the initial conditions of the costate 
variable is not a straightforward procedure. In practice, a 
common technique is to assume a final value of the costate 
variables and run the system backward in time from the ori- 
gin, using the optimal control law, until the state trajec- 
tory intercepts the initial state in the state plane. The 
idea of running time backwards is most useful for second- 
order systems, because all the switching is determined by a 
switching curve. It is of decreasing usefulness as the ord- 
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er of the system increases, as this would require switching 
surfaces in the phase space. Clearly, the ability of using 
second-order system techniques in the analysis of a high- 
order system is another advantage of the IMSC method. 

The differential equations for the elastic modes have 
the decoupled form 

5 r (t> + aj S r (t) = N r (t) r=l,2,...,n e (5.66) 

where n e is the number of elastic modes in the control 
scheme. Using the same state variable definition as in 
Eq.(5.6d), the modal state vector x r (t) satisfies the matrix 
differential equation 


jV^I 

I = 

r 

f— 1 

0 

1 

J 

| | 

f 0 ] 

I- + i 


( 

I ; 

!-w 2 o 

L r 


l X 2r ( C) J 

1 1 1 


N r (t) 


r = 1,2 


* n e (5.67) 


Here again, it is the control task to drive the modal plant 
form an arbitrary initial state given by 

x r C°) = [x lr (0) x 2r (0) 1 = U or l 0r l (5.68) 

to the state origin in the final time t^, or 
x r (tf ) =0 (5.69) 


while minimizing the modal fuel expenditure 
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f tf 

J (N ) = | N (t) | dt 

r r '0 


r - 1 , 2 , . . . ,n 


The control signal N r (t) is bounded in magnitude by 


(5.70) 


N (t) < N 

r r 


r 1,2,..., n 


(5.71) 


while the control bounds are defined by Eq.(5.23) 


The modal Hamiltonian associated with Eqs.(5.67) and 
(5.70) has the form 


f r ( ?r*?r’ N r’ C) " l N r (t) l + P lr (t) X 2r ( ° “ P 2r (t) “I X lr (t) 


+ P 2r (c) N r (c) 


r — 1,2,..., n 


(5.72) 


where the costate variables p^ r (i— 1,2) are the solution of 
the adjoint system differential equations 


lp, r <t)j L-l 0 . lp 2r (t> 


r = 1,2, . . . ,n 


(5.73) 


The minimization of the Hamiltonian with respect to 

N (t), leads to a control law similar to that obtained pre- 
r' 

viously for the rigid body mode solution, namely. 


N^Ct) = - N r dez (p 2r (t)} 


r = 1,2, ... ,n 


(5.74) 
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However, here, p 2r (t) is the solution to the costate 
Eqs.(5.73). Rewriting these equations in scalar form, we 
obtain 

P 2 r (t) + w 2 P 2 r (t) = 0 

P lr (t) = " p 2r (t) r= 1,2,... ,n g 

The solution to Eqs. (5.75a) becomes 

P 2r (t) = P r Sin ( V + 8 r ) r = 

where P r and 0 r are unspecified modal constants, 
boundary conditions were imposed on the costate vector 
p r (t). To determine a control law in Eqs. (5.74) that will 
drive the state to the origin, these constants must be spe- 
cified. 

To obtain the optimal solution to the state equations, 
we can rewrite Eqs. (5.67) in scalar form by substituting the 
optimal control law of Eqs. (5.74), so that 

X lr (t) + x ir (t) = N r (t) (5.77a) 

x 2r (t) = x ir (t) r=l,2,...,n e (5.77b) 

Equations (5.77) are subject to the boundary conditions 
(5.68). Considering the elements of control N r (t)=0 and 
N r (t)=±N r , one obtains the solutions 


(5.75a) 

(5.75b) 

(5.76) 
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\ 


x lr (t) = ^ sin o> r t + ? Qr cos u) r t 
r 


(5.78a) 


x 2r (t) «0r 


r r 

for N r (t)=0 and 


cos u t - ; sin oj t 
r or r 


r — 1,2,..., n 


(5.78b) 


’Or 


x, (t) = sin a) t + 

1 r oj r 


A ' 

1^0 r 


cos O) t + — 
r oj‘ 


(5.79a) 


X 2 r^ c) . for 

U) U) 


COS (i) t - 5 T 

- 1 or 


sin oj t r ~ 1 j 2 j ■ . i 


(5.79b) 


for N r (t)=±N r =A r 

State plane trajectories are obtained by eliminating the 

time t from Eqs.(5.78) and (5.79), with the result 

2 

r = 1,2, ...,n (5.80) 


x + 
lr 


fx 1 

2 

[L 

2r 

_ c2 + 

Or 

OJ 

— C 

s 0r 

; ui 

r J 


l r J 


( , > 

2 

r } 

2 

( K 1 

2 

(1 

A 

r 

X - — 7* 

lr oj z 
r 

V. 7 

+ 

X 2r 

OJ 

l r J 

= 

A 

5 0r 

r 

s* J 

+ 

5 0r 

OJ 

r 


for N r ( t ) =0 and N r (t)=A r , respectively. 


(5.81) 


Depending upon the modal initial conditions, Eqs.(5.80) and 
(5.81) represent families of circles in the modal state 
plane with coordinates x lr and x 2 r/w • The centers of these 

/H * 

circles lie at 0, +N /u|, -N r /w| on the x lr axis for the 

A A 

control actions 0, +N r , -N r respectively, as given by 
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Eqs.(5.74) and (5.76). We designate these families accord- 
ing to the related controls, such as the 0, the +N r and the 
-N family. Figure 5.7 shows these circles in the state 
plane, as well as an arbitrary state trajectory formed by 
the control sequence -N r for t 0 <t^t 1/ 0 for t!<t£t 2 , +N r for 
t 2 ,t<t 3 , 0 for t 3 <t£t 4 and +N for t 4 <t£t P . Clearly, any 
arbitrary trajectory in the state plane that terminates at 
the origin, must use circle X+ or X~ as the last portion of 
the trajectory. 

Based on this observation and utilizing the harmonic 
nature of the costate functions in Eqs.(5.73), one can con- 
struct qualitatively a trajectory pattern in the modal state 
plane. Figure 5.8 provides a description of the costate 
function for a certain control sequence, as well as the as- 
sociated trajectory. 

To recover the time f rdtn“the state plane, we use the defini- 
tion of the normalized time scale x=a)t, where t is the real 
time. It is shown in Appen. D that 

At=AB (5.82) 

where A& is an angle segment of a circular arc As along the 
traj ectory . 
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Figure 5.7: Modal State Trajectories Families for EM 



switch-off line I switch-on line 


Figure 5.8: Modal State Trajectory vs Modal Costate Function for EM 

(a) Modal Costate Function p, r (t) 

'(b) The Associated Modal State Plane Trajectories 
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Because the optimal controls, Eqs. (5.74) are functions of 
the costate functions p 2r (t), to construct the optimal state 
trajectory we utilize geometric information available from 
the description of the costate function. The fact that the 
trajectory is a path constructed from circular arcs given by 
Eqs . ( 5 . 80 ) and (5.81) is utilized as well. As can be seen 
from Fig. 5.8a, all the control portions are of the same du- 
ration, i.e., t AB =t CD =t E0‘ This observation, using 
Eq.(5.82), implies that the angles of the control segments 
in the state plane must be of the same magnitude, i.e., 

angle (AO'B) = angle (CO h D) = angle (E0~ 0) = (5.83) 

etc. These observations play an imnportant role in the con- 
struction of the state plane trajectory. A simple geometric 
approach shows that for Eqs. (5.83) to hold true in the state 
plane, the control sector angle a must be symmetric about 
the x 2r axis. In addition, it can be shown that 

a r = 3 r (5.84) 

These observations enable us to calculate the costate cons- 
tants P r and 0 r . 
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In Appen. D, two cases corresponding to different types of 
initial condition points in the state plane are investigat- 
ed. In the first, the point lies outside the control sector 
( > a/2 ) , 


as shown in Fig. D.2. For this case the costate amplitude 
constant is obtained as 
.2 1 21 


P = 
r 


1 - 


r n ur 
0 r 


2n N 
c r 


(5.85) 


where 


f 


r o - 


^0 r 


1 » 2 , . . . , n 


(5.86) 


is the length of the initial condition position vector in 
the state plane and n is the number of full non-zero con- 
trol strokes. The costate phase constant is obtained as 


a = tt - ’h = i - arctan 
r -r 


’0 r 


oj <; „ 

r Or 


1>2 , . . . , n 


(5.87) 


where the argument in Eq.(5.87) depends on the quadrant of 
location of the initial condition point H ( 5 0r , 5 o r /w r ) • 
the second case, the initial condition point lies inside the 
control sector ( $ < a 1 /2 ) in the state plane, as shown in 
Fig. D.3. For this case, the costate amplitude constant is 
obtained as 


P = 
r 


V* 


- 2 “T r 0 sin# 


fN 


4(n 2 -n ) 
c c 


r — 1,2, .. . ,n 


(5.88) 
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r 0 has the expression of Eq.(5.86). 

The costate phase constant for this case has the form 


9 r = il - ’P r + 5 r 3 H - arctan 


s 0r 


oj £ 
r or 


+ 5 r r = 1,2, 


.,n e (S.89) 


where, from Appen. D 


6 = arcsin 

r 


N 

— 7* COSt}> 

r 



r > 

N 

2 N ] 

1 

2 


r 

Z* 

r 

- 2ro sint}> 


_ 

rj 

1 w r 

; 


- arcsxn \ 


N , 

r a 

— COS ■=- 
2 
r 


■ 

r \ 

N 

2 

N 

ro + 

r 


- 2r 0 sin<p 


l r . 


r 


r = 1,2, ... ,n 


e (5.90) 


is the phase constant correction factor. In these expres- 
sions, $ r =$ r -ir/2 , as can be seen in Fig. 5.8, D.2 and D.3. 


5.3.1 Fuel Consumed and Final Time Selection EM 
1. Initial condition point lies outside the control sector 
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In view of Fig. 5.8, the modal fuel consumed can be calcu- 
lated as follows: 

F = N n c (5.91) 

r r c ar 


In this expression, n c tar is the total time for non-zero 
control applications. t ar is obtained by using the concept 
of normalized time, or 

t = ^ (5-92) 

ar u 

r 

Using Eg. (D2 . 6 ) , 


F 

r 


2n N 
c r 

arcsm 

to 

r 


ro u* 

2a N 
c r 


so that letting 


(5.93) 


ro 

— = pi 0 < pi < 1 (5.94) 

2n N 
c r 

the normalized fuel is obtained from Eq.(5.93) in the form 


r 1 

= — arcsm di 

“ r Pi 

Using again Fig. 5.8, 
be calculated as 


(5.95) 

the normalized time for solution can 


t = 4> + (n^ — l)ir + arcsin pj (5.96) 

Clearly, from Eq.(5.96) t->-~ as n or in view of Eq.(5.94), 
as Pi-^-O. However, the fuel in Eq.(5.95) reaches its abso- 
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lute minimum at the same time and it is a simple task to 
show that 


F 

r mm 


lim F 

Pl-0 r 


ro « r 


(5.97) 


A similar result was obtained in the case of RBM, where in 
order to obtain minimum fuel solution, the time had to ap- 
proach infinity. Here again, a realization of such a solu- 
tion is made by a selection of finite time for solution, in 
a similar manner as to the RBM case. 

Figure 5.9 provides a description of the normalized fuel 
given by Eg. (5.95) and the normalized time of Eq.(5.96) for 
the case 0 = 1 1 / 2 , n c =k 1 /p 1 for various values of k t , where 
ro ^2 

ki = — - — (5.98) 

2N 

r 

is a modal parameter. For this case, Eq.(5.96) has the spe- 
cial form 


fk. 



+ arcsin 


Pi 


and, in addition, 

lim x = kj 7T 
Pl-1 


(5.99) 


(5.100) 


2. Initial condition point lies inside the control sector 


4>± o/2 
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Following a similar pattern, the fuel consumed is obtained 
as follows 


rn - 2 — j rn sinO 

U) 

r 


= — < (2n -1) arcsin TT~1 — 

r w c 4(n^-n 


(n^-n ) 
c c 


-? ~ 4 


(5.101) 


Letting 


4(n 2 - n ) —t 
c c w 2 
r 


-- - P | 


0 ^ p 2 ^ 1 


(5.102) 


defining the modal parameter 


2 , r 

tq - 2 — j ro smp 


(5.103) 


and inserting Eq. (5.103) into Eq. (5.102), we obtain 


1 + (p 2 + 4k 2 ) ^ 


(5.104) 


Assuming that the correction factor 6 is small compared to 
other terms, Eq. (5.101) can be written as 


N / oj 
r r 


i + (p| + 4k 2 )a 


- 1 arcsin p 2 + 


(5.105) 
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Equation (5.105) represents an expression for the normalized 
modal fuel. 

In view of Fig. D.3a and using the same assumption for 6, 
the normalized time for solution is obtained in the form 

T = (m~n 2 ) + (n -l)a' + (n -l)(ir-a') 
c c 

“ N 

r 0 “ 2 ~7 r 0 sin* 

0J 

(5.106) 

Using the notation of Eq. (5.102), Eq.(5.106) can be reduced 

to 

r = * + (n c -l)ir 4- arcsin p 2 (5.107) 


= * + (n -1 )tt + arcsin 
c 


4(n- - n ) — - 
c C 1 0) 


Clearly, from Eq. (5.107) t-*-~ as n c -*«, or in view of 

Eq. ( 5 . 102 ) , as p 2 -»-0. However, the fuel in Eq. (5.105) reach- 
es then its absolute minimum 


F r min 5 U " F r ' IT (1 + 2k l> + * (5 - 108) 

p 2 ->-0 r 

Using the notation of Eqs. (5.102) and (5.103), the normal- 


ized time in Eq.(5.107) can be written as 


f 1 + (p“ + 4k 2 )‘ 


- 1 77 + arcsin P 2 
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Figure 5.10 provides a description of the normalized fuel of 
Eg. (5.105) and the normalized time of Eq.(5.109) for the 
case ^=0. For this case, Eq. (5.109) has the special form 
[l + (pf + 4k 2 )^ 


2p 2 

In addition, 


- 1 


Tf + arcsin p 2 


(5.110) 


lim t : y (1 + 4k?) 2 (5.111) 
P2-1 

After the time for solution has been selected then the modal 
fuel is calculated according to the process described in 
Figs. (5.9) and (5.10). 

The total modal fuel for the EM' s is obtained by summing all 
the modal fuels in Eqs.(5.95) and (5.105), depending on the 
different modal initial conditions. A numerical example for 
the fuel minimization problem is presented in Ch. VI. 




Chapter VI 
NUMERICAL EXAMPLES 


The control concepts presented in the last three chap- 
ters are illustrated by two numerical examples. The equa- 
tions of motion for the structure were derived by the finite 
element method, as described in Chapter II. Considering one 
element per beam member and four members per cable branch, 
the order of the system displayed in Fig. 1.1 is n=270. For 
the present examples, a reduced order model of the struc- 
ture, such as the one displayed in Fig. 2.1, is used. The 
model has 45 degrees of freedom. 


The following structure parameters were considered: 


Beam member ( aluminum) 
length of member 
cross-sectional width 
cross-sectional thickness 
mass density 
modulus of elasticity 
Poisson ratio 
rectangular cross-section torsi 


1=12 inches 
a=2 inches 
b=0.25 inches 
p=9. 93x10’ 2 lb/in 3 
E=10xl0 6 psi 

v=0 . 3 

factor k=0 . 3 
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Cable member ( steel ) 
length of member 
equivalent diameter 
mass density 
modulus of elasticity 
Poisson ratio 


1=20 inches 
d=0.0442 inches 
p=19 . 75x10' 2 lb/in 3 
E=24xl0 6 psi 
v=0 . 3 


Other structure configuration parameters, are as in Fig. 2.1 

An eigenvalue problem was solved for the reduced-order 
system. The first twelve natural frequencies are given in 
Tab. 6.1 and the associated modes are displayed in Fig. 6.1. 
The system is excited by initial velocities, generated by an 
impulsive force, as shown in Fig. 6.2. Mode participation 
factors for the first twelve modes are displayed in 
Fig. 6.3, in view of which we retain the first six modes for 
control, because no further benefits are gained by controll- 
ing higher modes. 

We propose to control the structure by the IMSC method, 
which requires as many actuators as the number' of controlled 
modes. Placement of these actuators is based on results ob- 
tained in Chapter IV, as well as on the location of the no- 
dal lines in the higher uncontrolled modes. Examining the 
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first six modes in Fig. 6.1, we observe that four of them 
have symmetric mode shapes, while the remaining two, have 
antisymmetric mode shapes. Actuators placed as shown in 
Fig. 6.4, meet the requirement of Eq.(4.12) for system cont- 
rollability. Moreover, actuators in this configuration 
barely excite the seventh, eight and ninth modes, because 
their locations are on, or close to nodal lines of these 
higher modes. Clearly, the given impulse at node 8 excites 
only the symmetric modes. The given actuators configuration 
represents a symmetric set. To this end, the actuators at 
nodes 1 and 11 must have identical characteristics to those 
at nodes 5 and 15 respectively. 
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Figure 6.1: 


The Twelve Lowest Modes of Vibration 
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impulsive force 
15 

10 



Figure 6.2: System Initial Conditions 
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Figure 6.4: Actuators Placement Configuration 
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6.1 LINEAR OPTIMAL CONTROL WITH QUADRATIC PERFORMANCE 
INDEX 


Implementation of Linear Optimal Control policy, using 
the IMSC method is based on the derivations of Chapter III. 
The control effort weighting factors R r , are chosen so that 
the lowest controlled mode reach an approximate settling 
time of 3 seconds for 5% of the initial amplitude. To ob- 
tain the appropriate modal damping ratio of c=0.33, we use 
Eq.(3.18) or Fig. 3.1, and select p=4.00. In addition, a 
structural damping ratio of £=0.05 was considered. The res- 
ponse parameters of Eqs.(3.15) are displayed in Tab. 6.1. 

The response is presented in the form of the transverse mo- 
tion in four selected nodes. Figure 6.5 displays the actual 
response at nodes 1,3,8, and 14 for the uncontrolled as well 
as for the controlled cases. Figure 6.6 presents the modal 
forces while Fig. 6.7 presents actuators forces versus time. 
Finally, assuming that the control forces are proportional 
to the rate of flow of fuel. Fig. 6.8 provides the total 
consumed fuel versus time. 
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TABLE 6.1 

System Parameters in the Linear Optimal Control Framework 


Mode 

r 

Natural 

Frequency 

0) 

r 

Damping 

Ratio 

5 r 

Modified 

Frequency 

0J 

nr 

Damped 

Frequency 

(i), 

dr 

Phase 

Angle 

' r 

i 

2.995780 

0.379679 

3.167651 

2.930451 

0.389449 

2 

3.259080 

0.379679 

3.446057 

3.188010 

0.389449 

3 

4.161381 

0.379679 

4.400122 

4.070633 

0.389449 

4 

8.371135 

0.379679 

8.851396 

8.188587 

0.389449 

5 

10.537362 

0.379679 

11.141903 

10.307579 

0.389449 

6 

12.932981 

0.379679 

13.674960 

12.650956 

0.389449 


7 15.402695 

8 18.471893 

9 19.938034 

10 23.332260 

11 32.394440 


12 


38.127151 





uncontrolled 


controlled 



node 1 


8 q 3 (t) 



node 3 


Figure 6.5: Actual Response at Nodes 1,3,8 and 14 




uncontrolled 


controlled 



node 8 



node 14 


Figure 6.5: (cont.) 

















151 


6.2 MINIMUM-FUEL PROBLEM 

Implementation of Minimum-Fuel Control policy, using 
the IMSC method is based on the derivations of Chapter V. 
It has been shown that the modal minimum fuel control is 
nonlinear of the "bang-off-bang" type. Most of the control 
techniques for this type of problem were developed as 
state (phase)- plane concepts. Therefore it proves conve- 
nient to show how the control task is achieved for each mode 
in the control scheme, using state-plane trajectories. 

Table 6.2 provides the system parameters for three cases of 
control bounds. The modal state-plane trajectories for 
these cases are shown in Fig. 6.9. Figure 6.10 presents the 
modal on-off control forces, while Fig. 6.11 presents actua^ 
tors quantized forces versus time. Finally, assuming that 
the control forces are proportional . to the rate of flow of 
fuel. Fig. 6.12 gives the total consumed fuel versus time. 

A comparison between fuel consumption in the linear op- 
timal control design. Fig. 6.8, and the minimum fuel design. 
Fig. 6.12, reveals a saving of about 20% for the latter. 
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TABLE 6.2 

System Parameters in the Minimum-Fuel Control Framework 


Case 

Control 

Bound 

M 

Mode 

r 

Modal State-Plane 
Initial Conditions 

Number 
of Control 
Strokes 
n 

c 

Modal 
Final 
Time 
tf (sec) 

r 0 

ip (rad) 



1 

0.039649 

0.0 

3 

2.27 



2 

0.0 

0.0 

- 

- 



3 

0.074399 

0.0 

7 

4.71 

1 

100 








4 

0.033788 

0.0 

15 

5.34 



5 

0.0 

0.0 

- 

- 



6 

0.003797 

-3.14 

6 

1.03 



1 

0.039649 

0.0 

2 

1.20 



2 

0.0 

0.0 

- 

- 



3 

0.074399 

0.0 

5 

3.16 

2 

180 








4 

0.033788 

0.0 

9 

3.08 



5 

0.0 

0.0 

- 

- 



6 

0.003797 

-3.14 

4 

0.53 



1 

0.039649 

0.0 

2 

1.20 



2 

0.0 

0.0 

- 

- 



3 

0.074399 

0.0 

4 

2.44 

3 

200 








4 

0.033788 

0.0 

8 

2.71 



5 

0.0 

0.0 

- 

- 



6 

0.003797 

-3.14 

4 

0.53 




Figure 6.9: Modal State-Plane Trajectories (Cases 1,2,3) 
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Figure 6.9: (coat.) 
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Figure 6.12: Total Consumed Fuel (Case 1) 



Chapter VII 

CONCLUSIONS AND RECOMMENDATIONS 


Recent designs of large structures require distributed- 
parameter models described by partial differential equa- 
tions. Exact solutions to these equations are limited to 
relatively simple geometries. However, in most practical 
problems, closed-form solutions are not feasible and system 
discretization is unavoidable. The discretization process 
converts the distributed-parameter system equations into a 
coupled multivariable linear eqations of motion. Control of 
the high-order model is not feasible due to computation li- 
mitations and errors that arise in modeling the higher modes 
by the discretization process. The fundamental problem of 
active control of flexible systems is the control of a 
large-order system, using a much smaller-order controller. 
The particular number of modes chosen to represent the trun- 
cated system model is often difficult to justify. One can 
measure the significance of the modes by computing their 
participation factors in generating the actual response. 

The increased order of models representing modern 
structures and the associated control problems have provided 
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the motivation for developing methods permitting analysis of 
such systems, and solutions to the control problems with 
little computation effort. Such a method is the Independent 
Modal Space Control (IMSC). Using this method, the 2n-order 
coupled system equations of motiom is reduced to n indepen- 
dent second-order modal space systems. This approach allows 
complete flexibility as to which, mode to control, because 
the control law is designed in the modal space for each mode 
independently. The ability of using second-order system 
techniques in the analysis of a high-order system is a dis- 
tinct advantage typical of the IMSC method and not possible 
for any other method. 

In this study, IMSC method was employed in the design 
of a linear optimal control law, as well as in the design of 
a minimum-fuel control law for a two dimensional domain 
problem. Due to the use of the IMSC method, some response 
performance indices, like settling time, can simply be met. 
At the same time, no control forces were applied to the 
modes that were not excited. These tasks are very hard to 
achieve, if not impossible, by any other method. 

An important factor in any control scheme is the place- 
ment of the actuators. Actuators placement becomes even 
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more critical when two dimensional domains are considered in 
conjunction with the IMSC method, especially if the system 
possesses some symmetry properties. In such cases, actua- 
tors placement must be carefully made to prevent the occu- 
rence of a singular mode participation matrix. Because ac- 
tuators location affect this matrix directly, knowledge of 
the modes shape permits a correct actuators placement. Re- 
lating to the modes shape of such systems, a method was de- 
veloped in this study as to how to recognize such singular 
cases just by examining the location, as well as the type of 
actuators used in the control framework. This method is 
very helpful in the design and placement of actuators, as 
well as of sensors when the IMSC method is applied to two- 
dimensional domain problems. 

The minimum- fuel problem is a very important one in 
space vehicles, especially in those designed for lengthy 
missions. The amount of fuel available for control is at 
times limited to such an extent that fuel economy is the 
predominant factor in control system design. The explicit 
dependence of the control on the instantaneous values of the 
state variables in a coupled high-order system, is a very 
complex task, and for most practical cases it remains an un- 
solved problem. In these cases, a trial and error process 
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is unavoidable for the state determination. Such a process 
may be reasonable for a low-order system (fourth order at 
most), but certainly not for a high-order system, where the 
computational difficulties are insurmontable . The control 
task is made considerably simpler by using the IMSC method, 
which obviates the complexity inherent in a high-order sys- 
tem. In this study, IMSC method was employed, and the dif- 
ficult minimum-fuel problem for a 2n-order coupled system 
was reduced to n fuel minimization problems for a decoupled 
set of n second-order modal space systems. In using the 
IMSC method, one must minimize modal cost functions instead 
of a global cost function. This decomposition allows the 
use of second-order system techniques in the control scheme 
rather than the multivariable system control techniques, 
which are hard to implement. 

In the derivation of the minimum-fuel control problem 
via the IMSC method, the condition number of the mode parti- 
cipation matrix plays an important role in determining actu- 
ators efficiency. Observing the characteristics of the ac- 
tuators quantised forces, one can easily determine this 
effectiveness to be less than 70%. It is suggested that 
further research be carried out on the dependence of the 
condition number on actuators placement and how to decrease 
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this condition number in order to increase actuators effec- 
tiveness. 

The fuel minimization technique depends on correct 
switching times in the process of driving the modes to the 
system origin. Variations in these times create an offset 
in reaching the final state. Further study of the sensitiv- 
ity of these switching times to parameters variation and 
their influence on the final state is strongly recommended. 
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Appendix A 

LAPLACE EXPANSION FOR DETERMINANTS 
See also Karlin [28], Chapter 0. 


The classical Laplace expansion for determinants has 
the form 


det A 


1 2 . . . n 

1 2 . . . n 


l (-D 


k=*l 


Uk + Jk> 


det A 


1} ±2 • • • i 

Ul J2 ... j p J 


fii i4 ...i 


• det A 


n-p 




(A.l) 


where 


ii it ...i, 


\2 1 3z 


Pj 


denotes the minor of the matrix A obtained by the selected 
rows i l ,i 2 ,...,i and selected columns ji/jz, • 


In the sum, there are ( n ) terms as the number of selec- 

P 

tions of ordered p indices. The sum is extended over all 

p- tuples l£j i <j z < . . • <i Sn, where j l , j 2 , . . . , j ' denote the 

p v 
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complementary set of indices to j 1 , j 2 , . . . , j p in l,2,...,n, 
arranged in natural order, and similarly i'j , i' z , . . . , i^_ p de- 
note the complementary set of indices to i 1 ,i 2 ,...,ip in 
1 , 2 , . . . , n. 

Example 1 

Expand the determinant 


det A = 


ail ai2 
a 21 a 22 
a 3 1 a 32 

a 4 1 a 42 


a l 3 

a l 4 

a 23 

a 24 

a 3 3 

a 34 

a4 3 

a 44 


(A. 2) 


using Laplace expansion. 

For the expansion, we will select two fixed rows, for exam- 
ple, the second and the third rows. The number of terms in 
the expansion will be 

< p > = < 2>= 6 

Using Eq. ( A. 1 ) 
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det A = + 


+ 


a 2 1 

a 22 

a 3 1 

a 32 

a 2 1 

a 24 

a 31 

a 34 

a 22 

a 24 

a 3 2 

a 34 


a l 3 
a 4 3 
a l 2 
a 4 2 
a U 
a 4 1 


a l 4 
a 44 
a l 3 
a 4 3 
a l 3 
a 4 3 


a 2 1 
a 3 1 
a 22 
a 32 
a 2 3 
a 3 3 


a 23 
a 33 
a 23 
a 3 3 
a 24 
a 34 


a l 2 
a 42 
a ll 
a 4 1 
a ll 

a 4 1 


a 14 
a 44 
a l 4 
a 44 
a l 2 
a 4 2 


Example 2 


(A. 3) 


We attempt to control the first six modes in a n=45 
order system, where the actuators are placed as shown in 
Fig. A.l. 


The mode participation matrix U C B for this case is 


mode 


shape 

4- 

17 

18 

23 

24 

29 

30 

S 

a l 

a 2 

a 3 

0 

a l 

-a 2 

A 

bi 

b 2 

0 

b 4 

-bi 

b 2 

S 

a i 

c 2 

c 3 

0 

ci 

-c 2 

S 

dl 

d 2 

d 3 

0 

dl 

-d 2 

A 

e l 

e 2 

0 

e 4 

” e l 

e 2 

S 


h 

f 3 

0 

fl 

"f 2 . 


■*- actuator # 



(A. 4) 




Figure A.l: Beam Lattice with Symmetric Actuators Configuration 
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where "S" and "A" denote symmetric and antisymmetric modes, 
respectively. Here, p+h=2+l=3 and s=4. Hence, p+h<s and 
the proposition states clearly that for this case, the mode 
participation matrix is singular. It is an easy task to 
show this, using rows and columns elementary transformation. 
We change the rows order having first the symmetric modes 
and then the antisymmetric modes (not necessary), so that 


s 

r ai 

32 

33 

0 

31 

-32 

s 

Cl 

c 2 

c 3 

0 

C 1 

-c 2 

s 

<*1 

d 2 

d 3 

0 

dl 

-d 2 

s 

*1 

f 2 

h 

0 

fl 

-f 2 

A 

*>1 

b 2 

0 

bu 

-bi 

b 2 

A 

l e l 

e 2 

0 

eu 

“ e l 

e 2. 


Substracting 1st column from the 5th and adding the 2nd to 


the 6th, 

we 

obtain 




f 3l 

a 2 

33 

0 

0 

0 ' 


C i 

c 2 

c 3 

0 

0 

0 


<*1 

d 2 

d 3 

0 

0 

0 


fl 

f 2 

*3 

0 

0 

0 


bl 

b 2 

0 

b 4 

-2b! 

2b 2 


L e l 

e 2 

0 

e 4 

- 2ei 

2e 2. 


Then, expanding the determinant of (A. 6), we obtain 
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dec ( ) 


which 


a l 

a 2 

a 3 

0 

0 

Cl 

c 2 

c 3 

0 

0 

di 

d 2 

d 3 

0 

0 

*1 

fz 

f 3 

0 

0 

e l 

e 2 

0 

~2ei 

-2e 2 


= 0 + 0 = 0 

T 

means that U C B is singular. 


a l 

a 2 

a 3 

0 

0 

ci 

c 2 

C3 

0 

0 

di 

d 2 

do 

0 

0 

fl 

f 2 

f 3 

0 

0 

b l 

b 2 

0 

-2b i 

-2b 2 


(A. 7) 



Appendix B 


MINIMUM FUEL CONTROL WITH UNSPECIFIED FINAL 
TIME: RIGID BODY MODE 


See also Athans [17]. 


To obtain the minimum- fuel solution and to indicate the 
set of initial states for which no minimum-fuel solution ex- 
ists, we divide the x x -x 2 state plane into four sets, using 
the 3f + and t~ curves defined by Eqs.(5.51) and (5.52) re- 
spectively, and the x lr axis. 


We recall the following equations to describe the state 

A \ 

plane trajectories for N r (t)=0, and N r ( t)=±N r =A r . 

Solution to the state equations 


x lr (t) 'hr* hr e 


(B.l) 


x 2r (c) ’ hr 


(B .2) 


for N ( t )=0 

T ' 


x ,r (t > 'hr + hr'+kh 


x (t) = ? + A t 

2r Or r 


for N r ( t)=±N r = A r 


( B . 3 ) 
(B.4) 
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The state plane trajectories are 

X 2r (t) = ^or = COnst » for N r (t) = 0 ( B . 5 ) 


X ir (t) = ^or 


_L_ r 2 

2A 


- + — x 2 

‘’Or 2 a 2r ’ 
r r 


for N (t) = ± N = A 
r r r 


( B . 6 ) 


The regions in the state plane are described in Fig. B.l. 
Now we show the regions of existence and non-existance for a 
minimum-fuel solution. 


B.l REGION R1 

If the initial state is in region R 1( given by 
Eq. (5.54a), or region R 3 , given by Eq. (5.54c) (note that the 
switching curve Z is not included), then a minimum- fuel so - 
lution for the control driving the plant to the origin does 
not exist. 


For any N r (t), we have from the second state equation 
^ 2 r ( t)=Nr( t) , which implies that 


x (t) 
2r 



N r (t)dx 


(Bl.l) 


from which 




Figure B.l: Control Dependence on Initial Condition Point for RBM 

Regions Ri , R 2 , R 3 , R 4 Are Defined by Eq. (5.54). 
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ft 



N ( T )d T 

= 

x„ (t) - £ 

't„ r 


2r ^Or 


(B1.2) 


Dividing the path to the origin ABCO into segments on which 
the control is a constant function of time, we can calculate 
the consumed fuel 


J (N ) = 
r r 


t' 


c o 

= |-6-| 


|N r (t)|dt + 


I N (t) |dt + 

I t 


r c f . , 

[N (C) | dt 

II r 


Or 1 


+ 1-6 -(-€) I + 0 - (-6) 


+26 
Or 1 


(B1.3) 


The elapsed time for this operation can be calculated as 
follows: 


on AB 
t’ - c ( 

on BC 


£ + 6 
or 

N 


t" - t' 


r + — — C 2 

Or 2 N r ^ 0r 


on CO 


C f - C ' 


N 


The final time tf is 


c + — =— 


^Qr + 6 . ~° r 

C f C 0 6 

N 


’Or 


(B1.4) 


(B1.5) 


(B1.6) 


(B1.7) 
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Clearly, the fuel in Eq.(B1.3) will reach its minimum in the 
limit as e-*-0. The required fuel, then, approaches the lower 
bound on fuel expenditure. However, as e ->0, the elapsed 
time in Eq.(B1.7) increases without bounds, tf-*-». Hence, 
the fuel required can be arbitrarily close to, but never at- 
tain the lower bound. This means that a minimum-fuel con - 
trol does not exist. 


B . 2 REGION R2 

If the initial state is in region R 4 given by 
Eq. ( 5 . 54d) , or region R 2 given by Eq. (5.54b), then the mini - 
mum-fuel solution for the control driving the plant to the 
origin is not unique . 

In this statement we excluded the case in which the in- 
itial state is on the switching curve % r . For such a case 
it is clear that a minimum-fuel solution exists and it is 
unique. [+N r } is applied to an initial state on as well 
as {-N r } to an initial state on Any other initial state 
can be represented by point D in Fig. B.l 
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Clearly, the fuel consumption would be the same for all 
paths DEO, DFGHO, DFGIJKO 

J r (N r ) = l^orl (B2.1) 

Hence, there are infinitely many ways of reaching the origin 
from D with the same amount of fuel. The associated control 
sequences with paths DEO, DFGHO and DFGIJKO are {0,+N r }, 
{0, +N r , 0, +N r } and [ 0 , +N r , 0, +N r , 0 , +N r } , respectively. The 
last two sequences, are certainly not from the set of candi- 
dates for the fuel-optimal control in Eq.(5.46), that minim- 
ize the Hamiltonian absolutely. In contrast to the same 
fuel consumption for those paths, the response time t f is 
different. The shortest response time among all minimum- 
fuel controls mentioned above is the one associated with 
path DEO and is given as 



Finally, the control sequence [-N r ,0,+N r ] along path DLMO 
cannot be fuel-optimal although it belongs to the set of 
candidates for the fuel-optimal control in Eq.(5.46). 

To show this, we assume that {-N r ,0,+N r } is fuel-optimal in 
driving the plant from its initial state at D to the origin. 
This control sequence is ploted in Fig. 5.3i, from which we 
can draw conclusions concerning the costate function P 2 r( t )- 
Using Eqs.(5.45), in view of Fig. 5.3i at t=0, we obtain 
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P 2 r (0) “ n 2 r > 1 P ir (0) = V (B2.3) 

Because N r (0)=-N r , substitution of Eq.(B2.3) into the Hamil- 
tonian Eq.(5.41) at t=0 yields 


H 


t=0 



lr or 


n (-N ) 
2 r r 


(B2.4) 


Because the final time tf is free and the Hamiltonian does 

not depend explicitly on time, the Hamiltonian must vanish 

on the extremal trajectory. The relation N r (l-n 2r )<0 im- 
• • • 
plies that m r Co r > 0* But, (5 o r ' S or ) 6 R 4 / so that £ or <0, from 

which we obtain n lr <0. In view of this result, we have 

ni r <0 and T) 2r >l such that from the costate equation, 

Eq. (5.45b), we obtain 

P 2 r (C) = n 2 r ‘ n ir C > n 2r > 1 for a11 c > 0 (B2.5) 


This corresponds to Fig. 5.3c with the associate control se- 
quence { -N r } , which contradicts the assumed sequence 
(-N r ,0,+N r |. Hence, the latter control sequence cannot be 
fuel-optimal . 



Appendix C 


MINIMUM-FUEL CONTROL WITH BOUNDED FINAL TIME: 
RIGID BODY MODE 


See also Athans [17]. 


C.l MINIMUM-TIME SOLUTION 


It will prove convenient to begin with the expression 
for the minimum-time solution. This expression appears fre- 
quently in the derivation for the minimum-fuel problem as it 
provides the lower bound for the final time. Figure C.l de- 
scribes a minimum- time optimal trajectory, where the optimal 

A A A /A 

control law is a sequence as [-N r ,+N r j or {+N r ,-N r j for an 
initial state in region R" or R + respectively. Regions R" 
and R + are the parts of the state plane created by the Y 
curve such that R" is to the right of this curve and R + is 
to the left. 


z l 


For an initial state in region R 

• i ~ 2 

£ + ± N t f 

Or Or 1 2 r 1 


Z') = £ - N ti 

2 or r 1 


0 = 2 ! + z 2 c 2 + y 


t2 


(Cl.l) 
(Cl. 2) 
(Cl. 3) 
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Figure C.l: State Plane Minimum-Time Trajectory for RBM 
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0 = 22 + t2 


(Cl. 4) 


where ti=t z -t B and t 2 =tf-t 2 . From Eqs.(C.l) (C.2) and (C.3) 
(C.4), we obtain 


t f - + ‘2 = t 0 + ^ + p£ + lf2l " 

N N 2 N 


(Cl. 5) 


For an initial state in region R + , we obtain in a similar 


t- - t 0 - —■ + 


2| 2 4£ ]* 

Or _ Or 

N 2 N 

r r 


{Cl. 6) 


C.2 SWITCHING TIMES TZ, TW 


Refering to Fig. 5.6, t is the first switch time and 
t is the second. The determination of these times is done 

W 

upon solving the following trajectories equation 


( Uor^or)eRiUR« ) 


path S-Z 


= 5 0r + V (t z - C 0> - 2 N r (C z " c <» : 


( C2 . 1 ) 


22 = « 0r - V C z - 


(C2.2) 


path Z-W 

W 1 =Zl+Z?(t - t ) 

1 1 - w z 


(C2.3) 
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w 2 = 2 2 


(C2.4) 


path W-0 


0 = Wl + w 2 (t f - t w ) + \ N r (T f - t w ) 2 


(C2.5) 


0 = w 2 + N (T c - t ) 
r £ w 


(C2.6) 


We express and t w in terms of the initial state (£or/Sor) 
and the prescribed final time Tf. Using Eqs.(C2.4) (C2.5) 


and (C2.6), we obtain 

2 


W 1 = 2 


i!i 

N 


From Eqs . ( C2 . 4 ) and ( C2 . 6 ) 


C w = T f + “ 
W N 

r 


and from Eq.(C2.2) 


(C2.7) 


(C2.8) 


t 2 - to ■ -r- <5 0r - z 2 ) 

‘ r 

Substituting Eq.(C2.9) into Eq.(C2.1), we obtain 

; 2 

f A Sr 22 

Z 1 = ^ 0r “ 

2N 2N 
r r 

From Eqs . ( C2 . 8 ) and (C2.9) 


2z 2 5 


t - c = — 
W 2 N 


Or 


N 


+ T f - t 0 


(C2.9) 


(C2.10) 


(C2.ll) 
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Substitution of Eqs.(C2.7) (C2.10) and (C2.ll) into 
Eq.(C2.3), yields a quadratic equation in Z 2 in terms of the 
initial state and the final time 


4 + [N r (T f - t 0 ) 



]z 2 + 




0 


(C2. 12) 


which has the solution 


z 2 " " 2 [M r (T ff - C 0> " ^0r ] 1 T\ [l V (T f " > “ *0^ 


- 4 


r 2 ‘ 

n e 

r or 2 


(C 2 . 13 ) 


In view of Eq.(C1.5), the terms in the curly brackets in 
Eq.(C2.13) add up to a positive value. Moreover, Eq.(C2.11) 

A • 

and knowledge that t w -t 2 >0 imply Z z >-l/2» (N r (Tf-t 0 )-$ 0 r ] / so 
that only the plus sign in Eq.(C2.13) is applicable. Sub- 
stitution of Eq.(C2.13) into Eq.(C2.9) and into Eq.(C2.8), 
yields 


1 

2 


' T f + 


“ 

T V ’ 0r 

r ; 

2 

i • 

o 

*1 

2 5 orl 

1*1 

Co + — - 
N 

r 

- 


N 

r 

j 

II 

j T 4 - 

• | 

’Or 

♦- 4. *• 4 » 

- 

r • \ 

T 

2 

.llor _ 

tc 2 1 


r 

c 0 + + 
N 

r 


1 f c 0 

N i 
r - 


N 

r 

N 2 
r - 

s 


(C2. 14). 


(C2.15) 


respectively. 
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C. 3 FIRST SWITCH TIMES LOCI 

Substitution of Eqs.(C2.7) and (C2.ll) into Eq.(C2.3), 
yields the equation for the loci of the first optimal switch 
time t z 

^Or _ (C3.1) 

“ — z 2 
N 

r J - 1 

( Zj , Z 2 )6 R 4 



C.4 BOUNDARIES FOR THE REACHABLE STATES 


Since t z and t w are real, the terms in the square 
brackets in Eq. (C2.14) or Eq.(C2.15) must add to be a posi- 
tive value. This condition yields the expression for the 
boundary to all modal initial conditions that can reach the 
origin in the prescribed final time bound, or 

< - -4“ [zf + 2d - t 0 ) N z 2 - (T - to) 2 (C4.1) 

4N 

r 

(Z x , Z 2 )e RxUR 4 Eq. (C4. l)can be modified to allow the determi- 
nation of the minimum final time bound for a given system, 
such that all the initial states in the control scheme will 
be located in the reachable zone. 


^ tg + max 

r-1,2 




(C4.2) 

l 0r )6R!URu 
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similarly. 


^ tQ + max 

r=l,2. 




- 4 


n 

r 1 J («or^Or )SR 2 TO 3 


(C4.3) 


Obviously, Eqs.(C4.2) and (C4.3) resemble the form of 
the minimum- time solution, in Eqs.(C1.5) and (Cl. 6), respec- 
tively. 


C . 5 PLOT OF THE MODAL STATE PLANE 


Figure C.2 presents a nondimensional plot of the modal 
state plane. The following scale factors have been used 


x, = x, N (T- - tn)' 
lr lr r f 0 


(C5 . 1) 


x 


2r 




to) 


(C5.2) 


Inserting Eqs.(C5.1) and (C5.2) into Eqs.(5.57) (C3.1) 


and (C4.1), we obtain the following nondimensional forms: 
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Figure C.2: Optimal Control Regions and Switching Curves in the 

Modal State Plane for RBM 
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1. For Eq.(5.57) - the curve that isolates modal initial 

states reaching the origin in less then the final time T 


or 


k l 2 + \ 

2 s or Or 


(C5.3) 


2. For Eq.(C3.1) - the family' of the first switching curves 
depending on the modal initial velocity 


| zl + <i - i or > z 2 


(C5.4) 


for £ or =const. 


3. For Eq.(C4.1) - the boundary curve for the reachable 

states 


£ 


or 




1 ) 


(C5.5) 


4. For Eq. (5.51) 
ing curve 


the curve, which is the final switch- 


(C5.6) 


Note that Eqs . (C5 . 3 ) - (C5 . 6 ) describe only one half of 
the plot in Fig. C.2, while the other part was obtained us- 
ing the symmetry about the origin. 
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Regions R',R°,R + in the plot designate the following 
control sequences to reach the origin. From region 

1. R' : { -N r , 0, +N r } 

2. R° : {0,+N r } fourth quadrant 

{0,-N r j second quadrant 

R + : {+N r ,0,-N r } 


3 . 



Appendix D 

MINIMUM- FUEL CONTROL: ELASTIC MODE 


See also Flugge-Lotz [23]. 


D.l TIME SCALING OF THE STATE PLANE TRAJECTORIES 


Figure D.l displays the construction for one short in- 
terval. Because time increases in clockwise direction in 
the state plane, the radius line R rotates in a clockwise 
direction. A positve increment As along the solution curve 
and a positive increment A(5 as well, are likewise taken in 
the clockwise direction of increasing time. From the geome- 
try of the figure, we obtain 


Axi x 2 
As u)R 


(Dl.l) 


Using the relations A(5=As/R, Ax!=x 2 At and the definition for 
normalized time t=ut, where t is real time, we obtain 


At * A8 


(D1.2) 


Equation (D1.2) implies that the normalized time is equal to 
the sector angle of the circle or circle approximation in- 
crement As of the state plane trajectory. 
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D. 2 DETERMINATI ON OF THE COSTATE CONSTANTS FOR INITIAL 
CONDITION POINT OUTSIDE THE CONTROL SECTOR 


The determination of the amplitude P and phase angle 0 
in the costate equation p 2 (t)=P«sin(ut+0) for a point out- 
side the control sector, is based on the following relations 
from Fig. D.2a 


0~ A 2 


In ' 

2 

N 

r 

r 

w ■> 


“ 2 *0 ^ c °s 

tt a 

[2 " 2 J 


(D2.1) 



fN 1 

2 

N 


> 

0- B 2 » OB 2 + 

r 


- 2 OB — 5- cos 

(t + 

a 

2 


Because 0~A=0'B, we obtain 
N 

Arg = 2 — r sin cf 
u an 2 


(D2.2) 


(D2.3) 


Note that Ar 0 = r 0 -OB is the decrease in the distance of the 
state from the origin of the state plane during one nonzero 
control stroke. Clearly, each such control action decreases 
this distance in equal steps. Defining the integer n c as 
the total number of control actions, we must have 


n c Ar 0 = r 0 


(D2 . 4) 


to reach the origin, where 


r 0 


.* 2 1 £ 
^0 


Co + -r 

u oj*- 


(D2.5) 


Substitution of Eq.(D2.3) into (D2.4), yields 



193 



Figure D,2: Initial Condition Point Outside the Control Sector 

for EM 

(a) State Plane Trajectory 

(b) Costate Function 
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a 

sm y ~ 


2n N 
c 


(D2.6) 


so that the only restriction on n c becomes 


> Integer 


(D2.7) 


The solution to the costate equation is given in Eq.(5.73) 


p 2 = P sin (wt + 9) 


(D2.8) 


Writing Eq.(D2.8) for point B in Fig. D.2b, we obtain 


1 = P sin £ + f 

4L 2 


from which, the amplitude constant P is obtained as 


P =• - 1 - 

ot 

cos 2 


.21 27~a 


2n c NJ 


The phase angle 0 is obtained directly from Fig. D.2b. 
Clearly, 


e + U- f + a + f - | 


from which 


9 + 4 + -y =■ tt 


and because 


♦ + \ = 


(D2.13) 
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we obtain 

0 = IT - ip = 77 


arctan 


fo 

^0 


(D2. 14) 


where the argument in Eq.(D2.14) depends on the quadrant of 
the initial condition point (5o/5o/«)* 


D. 3 DETERMINATION OF THE COSTATE CONSTANTS FOR INITIAL 
CONDITION POINT INSIDE THE CONTROL SECTOR 


The determination of the amplitude P and phase angle 0 
in the costate equation p 2 ( t )=P« sin( wt+0 ) for a point inside 
the control sector, is based on the following relations from 
Fig. D.3a 


"= 2 = r § + IfT 


0 = 


2 0 N fir , 

- 2r 0 ZZ C0S [2 ' * 


CfB 2 = OB 2 + 


Because 0'H=0'B, 


N ] 2 

N I 

tt a' 

—t* 

- 20B — r cos 

jj 1 

[2 + 2"j 


4 - OB* - 2 ^ 

u CO 


rg sinij) + OB sin y 
As in Eq. (D2 . 4) , 


(n c ~l)Ar 0 = OB 


(D3.1) 


(D3.2) 


(D3.3) 


(D3 .4) 
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where r^-l is the number of full nonzero control strokes and 
Ar 0 has the expression of Eq.(D2.3). Substitution of 
Eqs.(D2.3) and (D3.4) into Eq.(D3.3), yields 


sin 2 y 


2 „ N 

r 0 - 2 — r r 0 sin<p 


4 (n^ 


n c } 


0) 


(D3.5) 


Using the same argument as for Eqs.(D2.8) and (D2.9), we ob- 
tain the amplitude constant P for this case as 


P = 


N 


1 - 


— 2 r r q sind) 


4(n^ - n ) 
c c t or 


2 


(D3.6) 


The phase angle 0 is obtained directly from Fig. D.3b. 
Clearly, 


9 + (nx -n 2 )+j-|^ = Tr 


(D3.7) 


From Fig. D.3a we obtain the relations 


n l + Yl + (f ~ <t, J 


n2+Y2+f + f- = 7 T 


(D3.8) 


from which 


n l ~ 0 2 ~ Y2 " Yl + 


(D3.9) 
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Substitution of Eg.(D3.9) into Eq.(D3.7), yields 

9+ Y2“Yi + y+'J> = tt 


(D3. 10) 


or 

0 = TT - ’ 4 ) + 6 


(D3 . 11) 


where 6=Y 1 -X Z , is a correction factor to the phase angle 8. 

and are obtained from the triangular relations in 
Fig. D.3a, or 


sin y i 


sin y 2 



(D3.12) 


( D3 . 13 ) 


It is easy to show that 6—0 for the case where the initial 
condition point is outside the control sector (0>a/2). 
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